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PREFACE TO THE FIRST EDITION. 

OF all the celestial bodies whose motions have formed 
the subject of the investigations of astronomers, the Moon 
has always been regarded as that which presents the 
greatest difficulties, on account of the number of inequali- 
ties to which it is subject; but the frequent and important 
applications of the results render the Lunar Problem one 
of the highest interest, and we find that it has occupied 
the attention of the most celebrated astronomers from the 
earliest times. 

Newton's discovery of Universal Gravitation, suggested, 
it is supposed, by a rough consideration of the motions of 
the Moon, led him naturally to examine its application to 
a more severe explanation of her disturbances; and his 
Eleventh Section is the first attempt at a- theoretical 
investigation of the Lunar inequalities. The results he 
obtained were found to agree very nearly with those 
determined by observation, and afforded a remarkable 
confirmation of the truth of his great principle; but 
the geometrical methods which he had adopted seem 
inadequate to so complicated a theory, and recourse has 
been had to analysis for a complete ' determination of the 
disturbances, and for a knowledge of the true orbit. 



VI PKEFACE. 

The following pages will, it is hoped, form a proper 
introduction to more recondite works on the subject: the 
difficulties which a person entering upon this study is most 
likely to stumble at have been dwelt upqn'at considerable 
length, and though different methods of investigation have 
been employed by different astronomers, the difficulties met 
with are nearly the same, and the principle of successive 
approximation is common to all. In the present work, 
the approximation is carried to the second order of small 
quantities, and this, though far from giving accurate values, 
is amply sufficient for the elucidation of the method. 

The differences in the analytical solutions arise from 
the various ways in which the position of the moon may 
be indicated by altering the system of coordinates to which 
it is referred, or again, in the same system, by choosing 
different quantities for independent variables. 

D'Alembert and Clairaut chose for coordinates the pro- 
jection of the radius vector on the plane of the ecliptic and 
the longitude of this projection. To form the differential 
equations, the true longitude was taken for independent 
variable. 

To determine the latitude, they, by analogy to Newton's 
method, employed the differential variations of the motion 
of the node aad of the Inclination of the orbit. 

Laplace, Damoiseau, Plana, and also Herschel and Airy, 
in their more elementary works, have found it more con- 
venient to express the variations of the latitude directly, 
by an equation of the same form as that of the radius 
vector. 

Lubbock and Pontocoulant, taking the same coordinates 
of the Moon's position, make the time the independent 
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CHAPTER I. 

INTRODUCTION. 

BEFORE proceeding to the consideration of the moon's 
motion, it will be desirable to say a few words on the law 
of attractions, and on the peculiar circumstances which enable 
us to simplify the present investigation. 

1. The law of universal gravitation, as laid down by 
Newton, is that " Every particle in the universe attracts every 
other particle^ with a force varying directly as the mass of the 
attracting particle and inversely as the square of the distance 
between them" 

The truth of this law cannot be established by abstract 
reasoning ; but as it is found that the motions of the heavenly 
bodies, calculated on the assumption of its truth, agree more 
and more closely with the observed motions the more strictly 
our calculations are performed, we have every reason to 
consider the law as an established truth, and to attribute 
any slight discrepancy between the results of calculation and 
-observation to instrumental errors, to an incomplete analysis, 
or to our ignorance of the existence of same of the dis- 
turbing causes. 

Of the last source of error there is a remarkable ex- 
ample connected with the discoverv of the ulanet Neptune, 
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which became known to us, as one of the bodies of our 
system,* solely by means of the perturbations it produced in 
the calculated places of the planet Uranus. These pertur- 
bations were too great to be attributed wholly to errors of 
instruments or of calculation ; and therefore, either the law 
of universal gravitation ceased to hold for a body so remote 
as Uranus, or else some unknown cause was disturbing the 
path of the planet. The first supposition was too repugnant 
to an astronomer of the nineteenth century to be entertained 
until all others had failed, and the second supposition led 
to the detection of Neptune. The distinguished names of 
Adams and Le Verrier will be for ever connected with the 
history of this planet, and their solution of the difficult in- 
verse problem: "Given the perturbations caused by an 
unknown planet, determine, on the assumption of the truth 
of Newton's law, the orbit and position of this disturbing 
body," will always be considered one of the triumphs of 
Mathematical Science. 

Another remarkable instance of the vindication of Newton's 
law, where it seemed at first to be at fault, will be found 
mentioned in Chap. VI. in connection with the motion of 
the Lunar Apogee. The error was owing to an incomplete 
analysis; but, when the calculations were carried out more 
fully, the strongest confirmation of the law was afforded by 
the exact agreement of the results with observation. 

We shull proceed to apply the law of gravitation to in- 
vestigate the circumstances of the moon's motion ; and then 
see how it will enable us to assign her position at any time 
when sufficient data have been obtained by observation. 

2. The problem would be one of extreme, if not insur- 
mountable difficulty, if we had to take into account simul- 



* It had been scon by J)r. Lamont at Munich, one year before its being known 
to be a planet. '* Solar System, by J. R. Hind." 



IN PRODUCTION. 3 

taneonsly the actions of the earth, sun, planets, &c., on the 
moon ; but fortunately the planets are so small or so distant 
that their action may be neglected at any rate to the order 
of our present approximation and the attraction of the 
earth, on account of its proximity, is very much greater 
than the disturbing* action of the sun notwithstanding his 
enormous mass. V r e may therefore treat the question as 
that of one body, the moon, revolving about another body, 
the earth, and continually disturbed by a third body, the 
sun. This is the celebrated problem of Three Bodies. 

Stated in this general form its exact solution has hitherto 
defied the powers of the analyst. The disturbing action 
of the sun can be expressed without difficulty in terms of 
the masses and distances, but the integration of the diffe- 
rential equations cannot be effected. If however, instead 
of retaining the whole disturbing force of the sun, which, 
as we have said, is small compared with the earth's action, 
we expand its expression in a series and neglect all very 
small terms, it is then found possible to obtain a solution. 

This breaking up of the expression for the sun's action is 
obviously equivalent to a breaking up of the sun's force; 
and if we afterwards wish to take into account the very 
small forces which have been thus neglected, or the small 
disturbances due to the action of the planets, the following 
principle will shew that we may do so by considering their 
effects separately, and that the algebraical sum of all the 
disturbances so obtained will, to an order of approximation 
far beyond that which we contemplate, be the same as the 
disturbance due to their simultaneous action. We shall find 



* Since the sun attracts both the earth and moon, it is clear that its effects on 
the moon's motion relatively to the earth, or the disturbing force, will not be the 
same as the absolute force on either body. The absolute force of the sun on 
the moon is more than double that of the earth (Art. 99), but the relative force 
does not exceed T y of the earth's action (Art. 25). 
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an application of this principle in the investigation of the 
parallactic inequality. 

Principle of Superposition of Small Disturbances. 

3. Let a particle be moving under, the action of any 
number of independent forces, 
some of which are very small, 
and let A be the position of the 
particle at any instant. Let two of these small forces m^ m^ 
be omitted, and suppose the path of the particle under the 
action of the remaining forces to be AP in any given time. 

Let AL\ be the path which would have been described 
in the same time if m } also had acted; AP V differing very 
slightly from A P, the disturbance being PP^ 

Similarly, if 7/z 2 instead of n\ had acted, suppose AP^ to 
represent the disturbed path, PP 2 being the disturbance, 
(AP, AP t , AP t2 are not necessarily in the same plane, nor 
even plane curves). 

Lastly, let A Q be the actual path of the body when both 
m l and m z act. Join P^ Q. 

The two disturbances PP 2 and PjQ, being due to the 
action of the same small force m 2 on the slightly different 
paths AP and AP^ must be very nearly equal both in 
magnitude and in direction. What difference exists between 
them must be a very small fraction of either disturbance, 
and may be neglected compared with the original path. 
We may, in fact, expect that this difference which is the 
disturbance due to w t of a disturbance due to m 9 will be 
of an order compounded of the orders of the two. Thus, 
if the disturbance due to m t and tn z be respectively of the 
second and fourth orders, the difference between P^Q and 
PP 12 would probably be of the sixth order compared with 
A P. Therefore P 1 Q may be considered parallel and equal 
to PP t . 
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Hence the projection of the whole disturbance PQ on any 
straight line, being equal to the algebraical sum of the pro- 
jections of PP^ and JP t Q, will be equal to the algebraical sum 
of the projections of the separate disturbances PP^ PP a . 

Next, let there be three small disturbing forces m^ m^ w? a . 
We may consider the joint action of the two ?/?. 2 , m^ as one 
small disturbing force ; therefore, by what precedes, the total 
disturbance along any axis will be the sum of the separate 
disturbances of m^ and of the system ?/? 2 , m.^ ; but this last is 
the sum of the separate disturbances of ? 2 and wz a : therefore 
the whole disturbance equals the sum of the three separate 
disturbances. 

4. This reasoning can evidently be extended to any 
number of forces ; and if x, y, z be the coordinates of the 
disturbed particle, </> (a;, y, z) any function of #, y, s ; the 
disturbance produced in < (#, y } z} will be 



\ * * * *.*. f* \y e 

a?, y, z] - ^ &c -f -^ 8y + ^ &, omitting (Sx}\ &c., 

where 8^=8^ + 8^+...= sum of disturbances along axis of x 

due to separate forces, 
...= : ........................ along axis of y, 

...= ........................ along axisof 0; 



therefore 

W^C 

d<b ^ d<f> 



or total disturbance equals sum of separate disturbances, 
which establishes the principle. 

Since <f> (a?, #, s) may be the radius vector, the latitude 
or the longitude of the disturbed body, it follows that the 
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total disturbance in any of these elements is the sum of the 
partial disturbances. 

5. But we must still prove another proposition, without 
which the problem, though so far simplified, would scarcely 
be less complicated than in its most general form. 

Newton's law refers to particles, whereas the sun, earth, 
and moon are large nearly spherical bodies, and it becomes 
necessary to examine the mutual action of such bodies. 
Now, it happens that, with this law of force, the attraction 
of one sphere on another can be correctly obtained, and 
leaves the question in exactly the same state as if they were 
particles. (Princip. lib. I. prop. 75.) 



Attractions of Spherical Bodies. 

6. Let P be a particle situated at a distance OP a 
the centre of a sphere of uniform 
density p and radius c. The 
particle being without the sphere, 



from 



a> c. 




Let the whole sphere be divi 
ded into circular laminae by planes 
perpendicular to OP. Let SQ be 
one of these. PS=x, PQ = z, and thickness of lamina 

Next, let this lamina be divided into concentric 
Let RS-r be the radius of one of these rings and 
breadth, L RPO = ; 
therefore r = x tan 0, 



= $x. 
rings. 
Sr its 



The attraction of an element R of this ring on the particle 

. .,. , mass of clement , ^^ , . . , 

P will be --- P7?2~~~~~ a ^ on S "* anc * the resolved part 

^ . ! n ~ . . mass of element 

ol tins along PO will be -- ^ -^ - cosfl. 

x* sec'0 



ATTRACTIONS 



But the resultant attraction of the whole ring will clearly 
be the sum of the resolved parts along PO of the attractions 
of its constituent elements ; therefore, 

O c\ t\ 

attraction of rinff -= 4-*- ' cos STTD sin OSx . 80 : 
x secv r x 7 

rcos-i- 

thercfore, attraction of whole lamina SQ=27rp8xl siu0. 

Jo 



Again, = x 2 -f c 2 (a x)* = 2ax - (a 2 - c 2 ) ; 
therefore zSz = aSx, 

and attraction of lamina = 2?rp \- ---- , & ] : 

r V a 2a 2 J J 

/. attraction of whole sphere = 2?rp \ " --- ,- -- - , r 
1 r (2a 6a 4 2a j 

(from 5? = a c to z = 



M 



where -3/= mass of sphere = - ** . 

o 

Hence, the attraction of the whole sphere is precisely the 
same as if the whole mass were condensed into its centre. 

COE. 1. The attraction of a shell (radius c and thickness 
Sc) will be obtained from the preceding expression by diffe- 
rentiating it with respect to c, and is 

, 47rpc 2 . Sc mass of shell 

attraction of shell = ~ = , , 

or a ' 

the same as if the mass were collected at its centre. 
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COR. 2. Therefore, the attraction of a heterogeneous sphere 
on an external particle loill be the same as if the whole mass 
were condensed into its centre^ provided the density l>e the same 
at all points equally distant from the centre, for then the 
whole sphere may be considered as tne aggregate of an 
infinite number of uniform shells. 

7. Let us now consider the case of one sphere attract- 
ing another. Suppose P in the preceding article to be an 
elementary particle of a sphere 3/', whose centre 0' suppose 
at a distance a from 0. Then, since action and reaction 
are equal and opposite, P will attract the whole sphere M 
just as it would do a particle of mass M placed at 0. The 
same is true of all the elementary particles which compose 
the sphere M ', therefore the sphere M' will attract the sphere 
M as if the whole mass of the latter were condensed into its 
centre ; but the attraction of the sphere M' on a particle 
is the same as if the attracting sphere were condensed iato 
its centre 0' ; therefore, 

Two spheres attract one another as if the whole matter of 
each sphere were collected at its centre. 

8. This remarkable result, which, as may be shewn, holds 
only when the law of attraction is that of the inverse square 
of the distance, or that of the direct distance 3 or a com- 
bination of these by addition or subtraction, reduces the 
problem of the sun, earth, and moon to that of three par- 
ticles. The slight error due to the bodies not being perfect 
spheres is here neglected, being of an order higher than 
that to which we intend to carry the present investiga- 
tion : this error, however, though very small, is appreciable, 
and when a nearer approximation is required, it becomes 
necessary to have &egard to this circumstance. (See Ap- 
pendix, Art. 109). 



CHAPTER II. 

MOTION RELATIVE TO THE EARTH. 

9. When a number of particles are in motion under their 
mutual attractions or other forces, and the motion relatively 
to one of them is required, we must bring that one to rest 
and then keep it at rest, without altering the relative motions 
of the others with respect to it. 

Now, first, the chosen particle will be brought to rest 
by giving it at any instant a velocity equal and opposite to 
that which it has at that instant ; secondly, it will be kept 
at rest by applying to it accelerating forces equal and 
opposite to those which act upon it. 

Therefore give the same velocity and apply the same 
accelerating forces to all the bodies of the system, and t.'-eir 
absolute motions about the chosen body, which is now at rest, 
will be the same as their relative motions previously. 

Problem of Two Bodies. 

10. As the sun disturbs the moon's motion with respect 
to the earth, it is important to know what that motion would 
have been if no disturbance had existed, or generally : 

Two bodies attracting one another with forces varying 
directly as the mass and inversely as the square of the dis- 
tance, to determine the orbit of one relatively to the other. 

Let M) M' be the masses of the bodies, M ' being the body 
whose motion relatively to M is required, r the distance 



10 LUNAR THEORY. 

between them at any time , and the angle between r and 
some fixed prime radius. 

. M 
Tbe accelerating force of M on M' equals -5- acting towards 

M 1 ' 
MI wbilc that of M' on M equals - - in the opposite direction. 

Therefore, by the principle above stated, we must apply to 
both M and M' accelerating forces equal and opposite to this 
latter force, and M' will move about M fixed, the accelerating 

force on M' being = /^w y , if JJL = M 4- M ' and r = - . 

, T d*u & 

Hence, d^^ u== K^ 

f?f) 
where h = r* -7- = twice the area described in a unit of time, 

at 7 

and integrating 

w = ^{l-f ecos(0-a)J, 

e and a being constants to be determined by the circumstances 
of the motion at any given time. 

This is the equation to a conic section referred to its 

A* 
focus, the eccentricity being e, the semi-latus rectum , and 

the angle made by the apse line with the prime radius a. 

By observations made on the actual path of the moon, 
which cannot differ very widely from its undisturbed orbit, 
we infer that this latter would be an ellipse with an 
eccentricity about ^. 

In the same manner, if the sun and earth be the two 
bodies considered, the relative orbit would be an ellipse with 
an eccentricity about $. 

11. The angle a between the radius vector and the 
apse line is called the tru* anomaly. 
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If n is the angular velocity of a radius vector which, 
moving uniformly, would accomplish its revolution in the 
same time as the true one, both passing through the apse 
at the same instant; then 7^-f-a is called the mean 
anomaly where is a constant depending on the instant 
from which the time is reckoned^ its value being the angle 
between the prime ndius and the uniformly revolving one 
when = 0. 

Thus, let JtfT be the fixed line or prime radius, 

A the apse, 

M' the moving body 
at time , 

MIA the uniformly re- 
volving radius at same time, the 
direction of motion being repre- 
sented by the arrow. 

And, let Ml) be the position of Mp when t = 0, 

then TJ/Z) = e and is called the epoch,* 

DM/j, = nt, 

TMA OL longitude of the apse ; 
therefore, mean anomaly == A M/J, nt -f e - a, 
true anomaly = AMM' = T MM ' - TMA = - a. 

12. To express the mean anomaly in terms of the true in 
a series ascending according to the 2>owers of e, as far as c". 

2?r 2 area 

n . .. . = 2-7T ~ .-- 

periodic time h 

h 




* We shall sec (Art. 36) that the introductior of the epoch t is avoided 
in the Lunar Theory by a particular assumption ; but in the Planetary it forina 
one of the important elements of the orbit. 
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therefore h = na* (1 - e 2 )*, 

dt r* 2 tt~eY 1 



h h {1+6008(0-00}* 



= i (i _ : y) {i 2e cos (0 - a) + 3e" cos* (6 - a)} 

= - (1 - 26 cos (0 - a) + $e f cos2 (0 - a)} ; 

therefore 7i + e = 6 - 2e sin (6 - a) -f f <? 2 sin 2 (0 - a), 
or (w + e - a) = (0 - a) - 2e sin (6 - a) + f e 2 sin 2 (0 - a), 
the required relation. 

13. To express the true anomaly in terms of the mean to 
the same order of approximation. 

6 - a = nt + s - a + 2e sin (0 - a) - f e 2 sin 2 (0 - a),..(l) ; 
.'. 9 a = nt + s a first approximation. 
Substituting this in the first small terms of (1), we get 
a = fl 4. s a -t- 2e sin (rc 4- s a) ... a second approximation. 

Substitute the second approximation in that small term 
of (1 ) which is multiplied by e, and the first approximation 
in that multiplied by e*, the result will be correct to that 
term, and gives 

- a = nt + e a -f 2e sin {nt -f e - a + 2e sin (rc + e a)} 

- |e 2 sin2 (?i + 3 -a) 
= n -f s - a + 2e sin (?z + e a) 

H- 4e 2 cos(?z+s a) sin (w-f- s a) |e* sin2 (w-f e- a) 
= rctf 4- e a 4 2e sin (nt + s - a) + f c? sin 2 (w* + e - a), 
the required relation. 
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The development could be carried on by the same process 
to any power of e^ the coefficient of e* would be found 

|f sin 3 (nt -f - a) - -] sin (nt + a), 
but in what follows we shall not require anything beyond ^ 2 . 

ProMem of Three Bodies. 

14. In order to fix the position of the mooji with respect 
to the centre of the earth, which, by means of the process 
described in Art. (9), is supposed brought to rest, we must 
have some determinate invariable plane passing through the 
earth's centre to which the motion may be referred. 

If the sun and earth were the only bodies in the universe, 
then the plane of the elliptic orbit (in which, according 
to the last section, the motion would take place) would be 
fixed, and might be taken as the plane of reference ; but, 
as soon as we take into account the action of the moon 
and planets especially the moon the plane ceases to be 
fixed, and some other plane must be found not affected by 
these disturbances. 

Theory teaches us that such a plane exists,* but as its 
exact determination can only be the work of time, the 



* See Poinsot, u Tht'orie et dtttrmitiation de I'rr/tiateur tin systeme solaire" 
where he proves that .in invariable plane exists for tho solar system, that is, n 
plane whose position relatively to the fixed stars will always be the same, 
whatever changes the orbits of the planets may experience; but as its position 
depends on the moments of inertia of the sun, planets, and satellites, and there- 
fore on their internal conformation, it cannot be determined a priori, and ages 
must elapse before observation can furnish sufficient data for doing KO <l 
posteriori. 

This result Poinsot obtains on the supposition that the solar system is a 
free system ; but it is possible, as he furthermore remarks, nay probable, that the 
stars exert some action upon it, it follows that this invnr m ' f tbte plane may itself be 
variable, though the change must, according to our ideas of time and space, 
be indefinitely slow and small. 
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following theorem will supply us with a plane whose motion 
is extremely slow, and which may for a very long period, 
and to a degree of approximation far beyond that to which 
we shall carry our investigations, be considered as fixed and 
coinciding with its position at present. 

In what follows the action of the planets is too small 
to be taken into account. 

15. TJie centre of gravity of the earth and moon describes 
relatively to the sun an orbit very nearly in one plane and 
elliptic; the square of the ratio of the distances of the moon 
and sun from the earth being neglected.* 




Let Sj E, M be the centres of the sun, earth, and moon, 
G the centre of gravity of the last two. Now the motion 
ot Q is the same as if the whole mass E+ J/were collected 
there and acted on by forces equal and parallel to the moving 
forces which act on E and M. The whole force on O is 
therefore in the plane 8EM; join SG. 

Let L SGM~ w, and let m be the sun's absolute force. 

i in 
Moving force on E " in ES, 



f HJT wi* . ^-(3 

moving force on M= ^ m Mo* 



* This ratio is about ^-, and, as we shall see Art. (21), such a quantity 
we shall consider as of the 2nd*>rder of small quantities, and its square therefore 
of the 4th order. Our investigations are carried to the 2nd order only. 
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15 



These applied to G parallel to themselves are equi- 
valent to 





m'.E.GE 


m'. 


M.GM. 


r1itpftinn 


GM 




jinn 


SE a 
m'.E.SG 


m'. 


tiM" 

M.SG 




GS- 


by the 
, triangle 
of forces. 




SE 3 h 




SM' 









but 



E M M+E 



GM~ GE~ ME' 
whence E. GE= M. GM= (M+ E) 

Therefore substituting, and dividing by M + E, we get 
accelerating forces* 

, GM.GE 

m. 



anu 



m 



GE ( 1 1 \ j- - r, 

T (s& ~ -or) m directlon GM > 

E M_\ 

E 3 + 6'iWV ............... 



'.SG f E M 



M + E 

Again, SE' = SG* + GE* + 2SG . GE cos o> 
,A *OE QE*\. 







' \ 


1 ~ (j UUOU/ T y/va M 


therefore 


1&-B&' 


(' 


Q-TY- COSO)J 


omitting quantities 


similarly 


1 1 

SM* ~ SG* 


(l 


3 ^ Jf COSG)l 
-f COS ft) 1 


to be neglected. 








y 



* In strictness it would be necessary, since we have brought S to rest to 
apply to both M and t and therefore to G, accelerating forces equal and opposite 
to those which E and M themselves exert on S ; but the mass of S is so large 
compared with those of E and M, that we may safely neglect these forces in 
this approximate determination of the path of G, the error being of a still 

higher order than that introduced by the neglect of f - J J . 
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Therefore the accelerating force in the direction GM 
Zm'.GMME 

'--- cosw 



- 3 (accelerating force of sun on G) -^- . q ~~ cosw 



= 0, according to standard of approximation adopted. 
And the accelerating force in G8 

m'.SG (E + M 3(E.GE-M.GM) 



-B& SG< C08a) } 

= -^ to the same approximation. 

Hence, the force on G is a central force tending to 8 
and varying inversely as the square of the distance ; there- 
fore the orbit of G about 8 is very approximately an ellipse 
with 8 in the focus ; and the plane of this ellipse is, as far 
as our investigations are concerned, a fixed plane. 

This fixed plane is called the plane of the ecliptic, or 
simply the ecliptic. 

16. We shall now suppose a plane through the earth's 
centre parallel to the ecliptic; and this, when the earth is 
brought to rest, will give us the fixed plane of reference 
we require (14); the ecliptic then making small monthly 
oscillations from one side to the other of our fixed plane. 

17. Since G describes an ellipse relatively to the sun, 
the sun will describe the same ellipse relatively to G ; but, 
as seen from the earth, the orbit of the sun will be slightly 
different, as we shall now shew. 

First : The sun will have a latitude, that is, will be out 
of the plane of reference ; and this latitude will be of the 
same name (north or south) as that of the moon, and 
deducible from it when the ratio of the distances of the two 
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bodies from the earth, and that the masses of the earth 




and moon are known. For if S'EM' be the fixed plane 
through E, and #', G\ M' the projections of #, <7, M, then 
SES' is the sun's latitude, and MEM' is the moon's; arid 
the two bodies are obviously both on the same side of 
the plane. 

. , , . .. 88' GG' EG. sin (moon's lat.) 
sin (suns lat.) ==->^ == -77,7 = ----- 7 ^ t , --- - 
v ; ES ES Eti 

M EM . . , . . 



Now, it is known that M is about ^ of E, 



Al - . , , i x\ sin (moon's lat.) . 

therefore sin (sun's lat.) = -^ ---- ' nearly. 



And as the moon's latitude never exceeds 5 9', the sun's 
latitude will always be less than 1". 

Secondly, as to the sun's longitude: let EV be the 
direction of the first point of Aries, that is, a fixed line 
in our plane of reference from which the longitudes of the 
bodies are reckoned. Then TJ5/S" = 0' is the sun's longitude 
as seen from E- but, as seen from (7, the sun's longitude 
would be the angle between G 8 and a line through G (not 
drawn in the figure) parallel to ET. 

The difference between these two longitudes is the angle 



A 1 

a PP roxiraatel 7 5 
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therefore sinES'G' never exceeds ^{w an( * EffG 9 is a 
small angle not exceeding 7". 

Also US' - B 1 G' < EG 1 < <, S' G'. 

18. Now, by assuming the longitude and distance of the 
sun as seen from E to be the same as when seen from G, 
we commit the above small errors ii the position of S; 
that is, we assume the sun to be at 8" instead of $, 8' 8" 
being drawn equal and parallel to G'E. If our object were 
the determination of the sun's position, it would be neces- 
sary to take this into account ; but when investigating the 
motion of the moon, the disturbing action of the sun will 
not, on account of his great distance, be appreciably altered 
by supposing him to be at S" instead of P>, and the errors 
thus introduced are much too small to be considered. 

llence we may assume that the motion of the sun about 
the earth at rest is an ellipse having the earth for its focus, 
and its equation 

u' = a'{l+e' cos(0'-f)}, 

wh^.re is the longitude of the sun's perigee and e = ^ is 
the eccentricity, and we are safe that no appreciable error 
will ensue in the determination of the moon's place.* 

* That is, as far as the three bodies alone are concerned; but, since the 
attractions of the planets may, and in fact do, disturb the elliptic orbit of 
the sun about 6?, the same cause will disturb the assumed orbit about E. -A 
remarkable result of this disturbance is noticed in Appendix, Art. (108). 




CHAPTER III. 

RIGOROUS DIFFERENTIAL EQUATIONS OF THE MOON'S MOTION 
AND APPROXIMATE EXPRESSIONS OF THE FORCES. 

19. The earth having been reduced to rest by the process 
described in Art. (9), let us take its 2L 
centre E for origin. Let r, be the 
polar coordinates of the projection M' 
of the moon on the fixed plane of re- 
ference, being the longitude VEM' measured from a fixed 
line .Z?T in that plane. Let s be the tangent of the moon's 
latitude MEM ', so that MM' = rs. 

Next, let the accelerating forces which act Upon the moon 
be resolved into these three : 

P parallel to the projected radius M 1 E and towards the earth, 
T parallel to the fixed plane, perpendicular to P, and L* the 

direction of increasing^ 
8 perpendicular to the fixed plane and towards it. 

The equations of motion will be 



20. These three equations for determining the moon's 
motion take the time t for independent variable ; but it will 
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be more convenient in the following process to consider the 
longitude as such, and our next step will be to change the 
independent variable from t to 0. 
From (i) we get 



^^_ m 

dt dt( dtr dt> 

therefore ( r 2 -~Y= h* + 2 JTr'dff, 



= H 2 suppose, whence H -^ = 2V 3 . 



Therefore ^ = -- = Hu 2 . if u = - , 

dt r 7 r 1 

dt 1 1 / ^ r T , 



an important equation connecting the time t with the longi- 
tude 0. 

, . <& dr dO TT du 



ffr d du\ d du 



Sub&litute in (ii), 
therefore 



TJ. 
H 

jjj r>rp 

and substituting for H 2 and for H ^ their values l?-\ 2 ^ 

(its J U 

T 

and -5 , then dividing by AV, and transposing, 

d'u P T du d'u \ T 



du \ au/ 

..du TT dH 



this is called the differential equation of the moon's radius 
vector. 
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Lastly, 

\ ( ds du 1 

ds du\ 

di}- s Te) 

dH / ds du\ 

^--- - 

from(^l), 

.LI c r> n TT* 

therefore Ps- fi-JZ 

therefore 

d*s __ Ps-S T ds 

dP + *-'lfir-T?u*d0 
this is the differential equation of the moon's latitude. 

2 1 . If these three equations could be integrated under 
these general forms, the problem of the moon's motion would 
be completely solved ; for as only four variables w, 0, 5, and 
t are involved (the accelerating forces P, I 7 , and S beinr; 
functions of these), the values of three of them, as w, 0, 5, 
could be obtained in terms of the fourth t ; that is, the radius 
vector, longitude, and latitude would be known correspond- 
ing to a given time. 

But the integration has never yet been effected ; except 
for particular values of P, T, and S ; and the method which 
we are in consequence forced to adopt is that of successive 
approximation, by which the values of w, 0, and 5 are 
obtained in series in which the terms proceed according to 
ascending powers of certain small fractions. Some one 
fraction being chosen as a standard with which all other 
are compared, the order of the approximation is esteemed 
by the highest power of the small fractions retained. 
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It is usual to consider ^ as a small fraction of the first order, 
consequently -fa of ^V^Too > s ......... ............. second ...... 

BsW s ....................... third ...... 

and so on, other fractions being considered as of the 1st, 
2nd, &c. orders, according as they more nearly coincide with 



22. It is necessary therefore, before we can approximate 
at all, that we should have a previous knowledge (a rough 
one is sufficient) of the values of some of the quantities 
involved in our investigations; and for this knowledge we 
must have recourse to observation. 

We shall therefore assume as data the following results 
of observation : 

(1) The moon moves in longitude aoout thirteen times 
as fast as the sun. The ratio of the mean motions in longi- 
tude represented by m is therefore about T ^, and may be 
considered as of the 1st order.* 

(2) The sun's distance from the earth is about 400 times 
as great as the moon's distance. 

Hence the ratio of the mean distances = ^fa is of the 
second order.f 

(3) The eccentricity e' of the elliptic orbit which the 
sun approximately describes about the earth is about ^ ? 
and this, approaching nearer in value to ^V than to ? -J- , 
will be considered as of the 1st order. 

(4) During one revolution, the moon moves pretty ac- 
curately in a plane inclined to the plane of the ecliptic at 

* This approximate value of m is easily obtained j -the moon is found to 
perform the tour of the heavens, returning to the same position among the 
fixed stars, in about 27$ days; the sun takes 3f5J days to accomplish the 
same journey. 

f The distances of the luminaries may be calculated from their horizontal 
parallaxes, found by observations made at remote geographical stations. (See 
the Author's Astronomy, chap. XVI.). 
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an angle whose tangent is about -j^, and therefore of the 
1st order.* 

(5) Its orbit in this plane is very nearly an ellipse having 
the centre of the earth in its focus, and whose eccentricity is 
about equal to our standard of small fractions of the 1st order, 
^iz. "^V ? an( J tm " s w iH a ^ 30 b ver 7 nearly true of the projec- 
tion of the orbit on the plane of the ecliptic, f 

To calculate the values of the forces. 

23. We are now in possession of the data requisite for 
beginning our approximations, and we shall proceed to the 
determination of the values of 1\ T, and S in terms of the 
coordinates of the positions of the sun and moon. 

Let S, Ej M be the centres of the sun, earth, and moon, 




m, E, M their masses, 

E\ M'j the projections on the plane of the ecliptic, 

6r the centre of gravity of E arid M. 



* That the moon's orbit during one revolution is very nearly a plane inclined 
as we have stated, will be found by noting her position day after day among the 
fixed stars; and the sun's path having previously been ascertained in a similar 
way, the rules of Spherical Trigonometry will easily enable us to verify these 
statements. (A stronomy, p. 92) . 

f The elliptic nature and the value of the eccentricity of the moon's orbit may 
be found by daily observation of her parallax, whence her distance from the earth's 
centre may be determined : corresponding observations of her place in the heavens 
being taken, and corrected for parallax to reduce them to the earth's centre, will 
determine her angular motion. Lines proportional to the distances being then 
drawn from a point in the proper directions, the extremities mark out the form 
of the moon's orbit. 

A similar method applied to observations of the diameter of the sun will 
determine the eccentricity of its orbit. (J^/vwouy, p. 138). 
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The forces we Lave to take into account are, the forces 
which act directly on Jf, and forces equal and opposite to 
those which act on E] these last being applied to the whole 
system so that E may be a fixed point. 



Forces on M. 


Forces on E. 




E 


M 




Attraction of -E^iy ' m ME 


Attraction of-M= ,,^3. 


mEM 


M& 


MW 




attraction of S\m.MG . ,, 
' 1 W" m ^-^ 


attraction of S in. EG 
m' . "SB"*" 


in^Jf 


Wi . -rn \ ^"-*~ 

= ^i m MB 1 , ^ 








"~ SE' 2 >' 8H- 




|__ 11*11*^ tft fir /? 

equivalent to ' /SW 3 ^ c 


M-*-J lit . AJ \JT 

equivalent to j j&S' ' 


par 1 , to 08 



Therefore, the whole attraction upon J/, when E is 
brought to rest, is 

fMG EG\ 



and 



iriSG (~T - 



parallel to G8. 



These are the rigorous expressions of the accelerating 
forces on M, and they can be expressed in terms of the 
masses and coordinates of the bodies. 

Approximate values of P, T, 8. 

24. Our investigations will be carried only to the second 
order ; it will be sufficient therefore if, in the preceding, we 
neglect small quantities of the fourth and higher orders. 
Let /*-=.#+ Jf, 

MGM' = tan' 1 ^ = moon's latitude, 

E'T the direction of the first point of Aries, 



8G = / = -,; 



^ ff = longitude of sun, 



M ' E' = r = - ; WE' M' =0 = longitude of moon, 
/. SE'M' = - ff = difference of longitude of sun and moon. 
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Now, SM*=SM"+MM" 

= SGP+ GM" - 2SG.GM' coa SGM ' + MM'* 



therefore -g^ = 13 \ 1 + - cos (0 - 5') I ; 



for 5-5' or E'.M' differs from SGM' by less than 7', 

/ 6WI/V 
Art. (17), and f ,J is neglected, being of the fourth order. 



i i r ^ ^? T^' ^ 

Similarly, grr., = -3 j 1 - ,- cos(0 - fl')l ; therefore 
the accelerating forces on the moon are approximately 

+ ~ 8 (MG + G-E) ......... in direction ME, 



and ^-(GM'+ GE') cos(^- ff) ......... parallel to 

whence P= -= -+ ^' JKff\ cosMGM'- ^ M'E' cos'(^- 
T 1 r 



-/p (i - to - p B + i ^2 (^ - ff)} t 

= - % M'E' coa (6 - ff) sin (0 - ff) 



m'u' 3 s 



cos2 *-' 
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The differential equations in Art. (20), when these values 
of the forces are substituted in them, would contaiu the 
variables u and #'; but u is given in terms of 6' by the 
equation of Art. (18) u a' fl +e cos(0' f)}, and we shall 
find means to establish a connexion betVeen 6 and 0', which 
will enable us to eliminate the latter. 

25. Before proceeding further, it will be important to 
consider the order of the disturbing 'effect of the sun's action, 
compared with the direct action of the earth. If we examine 
the values of P, T, and S, it will be found that the most 
important of the terms containing m, which are clearly the 
disturbing forces since they depend upon the sun, are in- 

nv\ ** 

volved in the form ^ , while those independent of the sun's 
action enter in the form -5 . 



T 



We must therefore find the order of ^ compared with -^ , 



Now the orbits being nearly circular, and m the ratio 
of the mean motions, Art. (22), we have 

mean motion of sun periodic time of moon 
.mean motion of moon ~" periodic time of sun 

2-Trr* 



m 



therefore -^ : ~ :: m 3 : 1, 

or the disturbing force of the sun is of the second order 
when comnared with the direct action of the earth. 



CHAPTER IV. 

INTEGRATION OF THE DIFFERENTIAL EQUATIONS. 

SECTION L 

General process described. 

26. If in the differential equations we substitute tbe 
values of P, 2\ and 8 just obtained, we are led to the 
following forms: 



which cannot be integrated ; but it is found possible to obtain 
a solution which will be a first approximation, if we retain 
only the important terms on the right-hand side of each 
equation. This first approximation being substituted in the 
equations, the terms next in order of importance may be 
retained and a new solution becomes possible, which in its 
turn is the basis of a further approximation, and so on. 

First, neglect the disturbing force of the sun which is 
of the second order, and also the moon's latitude, which, 
as will be seen by referring to the expressions for the 
forces (Art. 24), will either enter to the Second power or else 
in combination with the disturbing force. 
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When this is done the equations become integrable, and 
values of u and s may be obtained in terms of correct to 
the first order; this value of u will then enable us to get 
the connexion between 6 and t to the same order. 

[Let us, however, bear in mind that the equations thus integrated are not 
the differential equations of the moon's motions, but only approximate forma 
of them; and it is, therefore, possible that the results obtained may be ap- 
proximations to the true solutions only for a limited range of values. 

Whether they are so or not, will appear by comparing them with what we 
already know of the motion from observation ; and this previous knowledge^ 
in the event of their not being approximations, will probably suggest such 
modifications of them as will render them so.] 

To proceed to the next approximation: We must sim- 
plify the right-hand members of our general equations by 
employing the values of u and s just obtained, retaining 
now terms of the second order which we had previously 
neglected. 

The integration can always be performed when the right- 
hand members of the equations are circular functions (sines 
or cosines) of or of its multiples ; and as our first approxi- 
mation will give us the values of u and s in that form, the 
new expressions for the forces will also be of the same 
character, and the equations be again integrable. 

Thus we shall obtain new values of w, *, t correct to the 
second order. These values, introduced in the same manner 
in the second members and terms of the next higher order 
retained, will lead to a third approximation, and so on, to 
any order; except that if we wish to carry it on beyond 
the third, the approximate values of the forces, given in 
Art. (24), would no longer be sufficiently exact, and \ve 
must obtain more correct values from the expressions of 
Art. (23). 

27. We must tare stop to notice a peculiarity in these 
equations, when solved by this process. We have said that 
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to obtain the values to any order, all terms up to that order 
must be retained in the second members ; but it may happen 
that a term of an order beyond that to which we are work- 
ing would, if retained, be so altered by the integration as 
to come within the proposed order. 

Such terms must therefore not be rejected, and we 
shall proceed to examine by what means they may be 
recognised. 

28. Suppose then that after an approximation to a certain 
order, the substitutions for the next steps have brought the 
equation in u to the form 



where the coefficient G is one order beyond that which wo 
intend to retain. The solution of this equation will be of 
the form 

w= ...... + G' cos(p0 + #)+ ...... 

G' being a constant which may be determined by substi- 
tuting this value of u in the differential equation, 

whence G' = - - & ; 

from which we learn that if p differs very little from 1, G' 
will be at least one order lower than #, and will come 
within our proposed approximation; and consequently the 
term G coa(p0 + H) must be retained in the differential 
equation. 

The equation of the moon's latitude being of the same 
form as that of the radius vector, the same remarks apply 
to it. 

29. If p is small, G' = G very nearly, and the term 
G cos(p0 -f- H] will not rise in importance in the value of u. 



30 LUNAR THEORY. 

We should therefore be led to reject it in the differential 
equation, except for the following reason : 

In finding the connexion between the longitude and the 
time (one of the principal objects of the Theory), we must 
use the equation 



Now, substituting for u its value in terms of 0, we shall 
find that such a term as G' cos (pO -f- H) in u gives rise to 
a term G" coz(p0 + II) in this equation, where G' and G' 
are of the same order, then 

~ = ...... G"Q 

G" 
whence = ...... Bin(j>0 + J7); 

G" 
therefore, when p is a small quantity of the first order, 

will be one order lower than G", and the term will have 
risen in importance by the integration. 

T 
But yet further, if such terms occur in cy- s , they will be 

tit U 

twice increased in value ; for they increase once in forming 

I T 

\ / r~3 ^5 an( i once again, as above, in finding t. 

30. We have, therefore, the following rule : 

In approximating to any given order to the values of u 
and s, we must, in their differential equations, retain periodical 
terms ONE ORDER beyond the proposed one^ when the coefficient 
of in their argument* is near unity. 

*The angle of a periodical term is its only variable part and is called 
*he argument. 
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In approximating to the value of #, we must, in the same 
equations, also retain terms ONE ORDER beyond the proposed 
one when the coefficient of the argument is near zero; and 

T 
when such terms occur in ;-, -, they must be retained TWO 

fi U 

ORDERS beyond the proposed approximation.* 

31. We shall here, for convenience of reference, bring 
together the equations and the approximate expressions for 
the forces : 

P T du /d' 



* Instead of the forces which really act on the moon, we originally substituted 
three equivalent ones. P, T t S j these again are, by the preceding expressions, 

P 
replaced by a set of others. For, we may conceive each of the terms in p- 3 , &c., 

to correspond to a force, a component of P, T 7 , or 8] each force having the 
same argument as the term to which it corresponds, and therefore going through 
its cycle of values in the same time. Now, by Art. (29), when the coefficient of 
6 in the argument is near unity, the term becomes important in the radius vector, 
and when near zero, in the longitude : hence, a force whose period is nearly tho 
same as that of the moon, produces important effects in the radius vector ; and a 
force whose period is very long will be important in its effect on the longitude. 
See Airy's Tracts, Planetary Theory, p. 78. 
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SECTION II. 
To solve the Equations to the first, order. 

32. Neglect all terms depending on the disturbing force, 
i.e. those which contain m '; such teriiis being of the second 
order. Art. (25). 

The latitude s of the moon can n n ver exceed the inclina- 
tion of the orbit to the ecliptic ; but this inclination is of the 
first order, Art. (22), therefore s is at least of the first order 
and s* may be neglected. 

wi. p T A A-S A 

Whence ~ 2 = ^; - = 0; - y- =0, 

and the differential equations become 



whence u = p {1 + e cos(0 a)} ; or, writing a for ~ , 

w = a {l+ecos(0-a)} 
and s==Asin(^~7) 

c, a, i, 7 being the four constants introduced by integration. 

33. These results are in perfect agreement with what 
rough observations had already taught us concerning the 
moon's motion Art. (22); for u-a [\ +ecos(0- a)} repre- 
sents motion in an ellipse about the earth as focus. 

And s = k sin (0 7) indicates motion in a plane inclined 
to the ecliptic at an angle tan^Ar. 

For, if TOM' be the ecliptic, 
M the moon's place, MM 9 an arc 
perpendicular to the ocliptic, then 
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0; and if TO be taken equal to 7, and OM joined 
by an arc of great circle, we have 

tan MM ' = ten MOM' sin OM' ; 
or s == tan MOM ' sin (6 - 7), 

which, compared with the equation above, shews that 



Therefore the moon moves in a plane passing through 
a fixed point and making a constant angle with the 
ecliptic. 

The values of a and 7 introduced in the above solutions 
are respectively the longtitude of the apse and of the node. 

34. What the equations cannot teach us, however, is 
the magnitude of the quantities e and k. For this we must 
have recourse to observation, and by referring to Art. (22), 
we see that e is about ^V and k about T 1 2 , that is, both 
quantities are of the first order. Their exact values cannot 
yet be obtained: the means of doing so from multiplied 
observations will be indicated further on. 

35. Lastly, to find the connexion between t and 0, the 
equation (<7) becomes, jT=0, 



_ __ _ 

dO ~ hu* ~ h(f [l + e cos(0 - a;}* " 

Now this is the very same equation that we Lad found 
connecting t and in the problem of two bodies. Art. (12), as 
we ought to expect, since we have neglected the sun's action. 
Therefore, if p be the moon's mean angular velocity, we 
should, following the same process as in the article referred 
to, arrive at the result 



which is correct only to the firs* order, since we have 

p 
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rejected some terms of the second order by neglecting the 
disturbing force. 

36. The arbitrary constant e, introduced in the process 
of integration, can be got rid of by A proper assumption: 
this assumption is, that the time t is reckoned from the 
instant when the mean value of 6 is sero.* 

For, since the mean value of 0, found by rejecting the 
periodical terms, is pt + &] if, when this vanishes, = (), we 
must have e = ; therefore 



correct to the first order.f 

37. We have now obtained three results, ZTj, 8^ 1? 
as solutions to the first order of our differential equations, 
and we must employ them to obtain the next approximate 
solutions; but, before they can be so employed, they will 
require to be slightly modified in such a manner, however, 
as not to interfere with their degree of approximation. 

The necessity for such a modification will appear from 
the following considerations : 



* When a function of a variable contains periodical terms which go through 
all their changes positive and negative as the variable increase? continuously, 
the mean value of the function is the part which is independent of the periodical 
terms. 

f We shall also employ this method of correcting the integral in our next 
approximation to the value of tf in terms of t ; and if we purposed to carry our 
approximations to a higher order than the second, we should still adopt the 
same value, that is, zero, for the arbitrary constant introduced by the integration. 
To shew the advantage of thus correcting with respect to mean values : suppose 
we reckoned the time from some definite value of 0, for instance when = 0; 
then, in the first approximation, 

= i + 2e sin (t - a) 

is the equation for determining the constant t t and in the second approximation, 
would be found from 

= + 2c sin ( - o) + \& sin 2 (t-a)+ , 

giving different values of at each successive approximation. 
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Suppose we proceed with the values already obtained; 
we have, by Art. (24), 



= a(t - JO- Ij {I + * cos(0-a)P 

= a-r ............. -f u4 cos(0- a)-f ...... ; 

and this being substituted in the differential equation (A) 
of Art (31;, gives 



_ ............ , 

the solution of which is 

u~a{l+e cos(0 - a)} + ...... + \AQ sinCfl - a). 

Our first approximate value u = a {1 -f c cos(0 - a)} is here 
corrected by a term which, on account of the factor 0, admits 
of indefinite increase, and thus becomes ultimately more 
important than that with which we started. Such a cor- 
rection is inadmissible, for our first value is confirmed by 
observation (22J to be very nearly the true one. The moon's 
distance, as determined by her parallax, i,s never much less 
than GO times the earth's radius; whereas this new value 
of ?/, when 6 is very great, would make the distance in- 
definitely small. On the same principle, we see that any 
solution, which comprises a term of the form ^0sin(0 a), 
cannot be an approximate solution except for a small range 
of values of ff. 

Such terms, * if they really had an existence in our gyf<tem, must end in its 

* destruction, or at least in the total subversion of itH present state ; but when 
'they do occur, they have their ongm, not in the nature of the differential 
' equations, but in the imperfection of our analysis, and in the inadequate repre- 
1 sentation of the perturbations, and are to be got rid of, or rather included in 

* more general expressions of a periodical nature, by a more refined investigation 
*than that which led us to them. The mature of tMa difficulty will be easily 

* understood from the following reasoning. Suppose that a term, such aa 
'a sin (AQ + B], should exist in the value of , in which A being extremely 
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1 minute, the period of the inequality denoted by' it would be *of great length ; 
' then, whatever might be the value of the coefficient a, the inequality would 
1 still be always confined within certain limits, and after many ages would return 
< to its former state. 

'Suppose now that our peculiar mode of arriving at the value of u led us 

* to this term, not in its real analytical form a pin f 10-f /?), but by the way of 

* its development in powers of 0, a -f fid 4- X0* -f &c. ; and that, not at once, but 

* piecemeal, as it were; a first approximation giving us only the term a, a 

* second adding the term flO, and soon. If we stopped here, it is obvious tha^ 
'we should mistake the nature of this inequality, and that a really periodical 
'function, from the effect of an imperfect approximation, would appear under 
'the form of one not periodical These terms in the value of t/, when they 

* occur, arc not superfluous; they are essential to its expression, but they lead 
'us to erroneous conclusions as to the stability of our system and the general 
' laws of its perturbations, unless we keep in view that they ore only parts of 
'series; the principal parts, it is true, when we confine ourselves to intervals 
' of moderate length, but which cease to be so after the lapse of very long times 

* the rest of the series acquiring ultimately the preponderance, and compensating 
'the want of periodicity of its riibt terms.-HiR JOHN HERSCHEL, Encyclopaedia 
Aktropolitana PHYSICAL ASTUOM>MY, p. 079. 

38. To cxtritcate ourselves from this difficulty, and 
to alter the solution so that none but periodical terms 
may be introduced, let us again observe that the equation 

- ,:r t + M= = a ? which gave the solution /j and thus led 

to the difficulty, is only an approximate form of the first 
order of the exact equation. Any value of w, therefore, 

which satisfies the approximate equation -^, + w = a to the 

first order, and which evades the difficulty mentioned above, 
may be taken as a solution to the same order of the exact 
equation. 

These conditions will be satisfied if we assume 

u = a {1-f e cos(c0 - a)} ; 
then ,-- . 4- u = a 4- ae(l - c 8 ) cos(e# a) 

uu 

=35 a to the*first order, 
provided 1 c 7 be of the first order at least. 
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39. The introduction of the factor c in the argument is 
an artifice due to Clairaut ; and its effect, as we shall see 
in Art (66) is equivalent to supposing that the apse of the 
moon's orbit is not fixed. Now, although the observations 
recorded in Art. (22) did not suggest this, we must bear in 
aiind that they were extremely rough, and carried on only 
for a short interval; but when they are made with a little 
more accuracy, and extended over several revolutions of the 
moon, it is soon found that both her apse and the plane 
of her orbit are in constant motion. 

Clairaut was fully aware of these motions, and there is 
no doubt that he was led to the above form of the value 
of u by that consideration, and by his acquaintance with the 
results of Newton's ninth section, which, when translated into 
analytical language, lead at once to the same form.* 

We might, therefore, taking for granted the results of 
observation, have commenced our approximation at this step, 
and have at once written down u = a {1 4- e eos(c# a)} ; but 
we should, in so doing, have merely postponed the difficulty 
to the next step, since there again, as we shall find, the 
differential equation is of the form 

-f u = a function of 0. 



* Newton has there shewn, that if the angular velocity of the orbit bo to 
it of the body as G - F to G, the ad 
the original force being /u" 2 . Therefore 



G'* 
that of the body as G Fto G, the additional centripetal force is " 



F 

where ^ is the same as our c. 
G 



u = a 1 4* e cos 



( - a H , 
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the correct integral of which would be 



and this would at the next operation bring in a term with 
6 for a coefficient, which we now know must not be. We 
shall, therefore, hereafter omit such terms as A cos (5 JB) 
altogether, and merely write 

- a;} -f ............. 



40. The value of c will be found more and more' cor- 
rectly at each successive approximation, by always assuming 
these terms as the first two terms of the value of w, then 
substituting in the differential equation and equating coeffi- 
cients. It will thus be found that ae(l c 1 ) must equal the 
coefficient of cos(c# a) in the differential equation; and this 
will enable us to determine c to the same degree of approxi- 
mation as that of the differential equation itself. See 
Arts. (52) and (103). So far, all that we know about c is that 
it differs from unity at most by a quantity of the first order. 

41. In carrying on the solution of s, the same difficulty 
arises as in w, and it will be found necessary to change it into 



g being a quantity which differs from unity at most by a 
quantity of the first order. 

The introduction of g is connected with the motion of the 
node in the same way as that of c is with the motion of the 
apse (85j ; and the value of g will be determined by a process 
similar to that explained above for c. Arts. (51) and (102 j. 

42. The connexion between 6 and t will also be modified 
by this change in the value of w, 

dt _ 1 __ 1 _ 
dO ~ ha 9 [1 + e cos(c0 - a;] 5 " 
d.ct 1 1 



or 



d.c0 ha* {1 + e cos(c0- a)}' * 
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Here c9 and ct hold the places which and t occupied 
in (35); therefore 

cO = cpt 4- 2e sin (cpt - a), 
or 6 = pt + 2e sin (cpt a), 

6 

+o the first order, since - = e to the first order. 
c 

43. Since the disturbing forces are to be taken into 
account in the next approximation, we shall have to- use 
the value of u found in (18), which is 

n' = a'{l+' cos(0'-f)}: 

but this introduces ff ; we must therefore further modify it 
by substituting for & its value in terms of 0, and it will 
be found sufficient, for the purpose of the present work, 
to obtain the connexion between- them to the first order, 
which may be done as follows : 

Let m be the ratio of the mean motions of the sun and 
moon, 

p', p their mean angular velocities; .\y> = mp, 
p't + ftipt ...... mean longitudes at time , /3 being the su n ' 

longitude when = 0, 
ff, ...... true longitudes at time f, 

f, a ...... longitude of perigees when t = ;. 

therefore & f = sun's true anomaly, 

and pt -f # - f = ...... mean anomaly. 

But, by Art. (13), 

true anomaly = mean anomaly + 2e sin (mean anomaly) + &c. ; 
therefore 0'= p't+0 +2e'sin(/ + /3-?) + ...... 



to the first order j 
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because pt~0 2e sin (c0 a) to the first order by 
Whence u = a {1 + e cos (inO -f /3 - f)} to the first order. 

44. The values of sin 2 (0- 0') and cos 2 (0-0') can also 
be readily obtained to the same order : 

sin 2 (0-0') 

= sin {(2 - 2m) 0-20- 4e' sin (mO + fc- )} 



Similarly, cos 2 (5 - 0') = cos {(2 - 2wt) - 2^8} 
- 2e f cos {(2 - m) - - ?} + 2e cos {(2 - 3.) - 3/3 + "}. 
The first term of each of these is all we shall require. 

SECTION III. 

To solve the equations to the Second Order. 

45. Let us recapitulate the results of the last approxi- 
mation. 

u = a {1 -f e cos(c0 - a)], 
u' = a {1 + e' cos(0 + /3 - f )}, 
7), 



These values must now be substituted in the expressions 
for the forces, retaining terms above the second order, when, 
according to the criterion of Art. (30), they promise to be- 
come of the second order after integrating. 
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The differential equations (A) and (B) will then assume 
the forms 



and the integration of these will enable us to obtain u and s 
to the second order; after which, equation (C) will give the 
connexion between 6 and t to the same order. 

46. The quantity -,-^ r , which we shall meet with as a 

J ha 7 

coefficient of the terms due to the disturbing force, can be 
replaced by m' 2 a, m being the ratio of the mean motions of 
the sun and moon. 

So long as we neglected the disturbing force, h and a 
had determinate values: they belonged to the ellipse which 
formed our first imperfect solution, and would therefore bo 
known from the circumstances of motion in that ellipse at 
any instant ; h being double the area described in a unit 
of time, and a the reciprocal of the semi-latus rectum. It 
would consequently be impossible to assume any arbitrary 
connexion between them. But, when we proceed to a second 
approximation and introduce the disturbing force, there is 
no longer a determinate ellipse to which the h and a apply : 
the equation /<& = tfa of Art. (32) merely shews that a and h 
must refer to some one of the instantaneous ellipses which 
the moon could describe about the earth if the disturbance 
were to cease, and we are at liberty to select any one of 
these which will allow us to proceed with our approximation. 
The particular ellipse is determined by the above assumed 






relation -^5-5 = w 2 a, and the selection is suggested and 
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justified by the following reasoning : 

_JP' _ av ^rage period of moon about earth 
~~ P ~~ period of sun about earth * 

but, since the instantaneous ellipses are nearly circles, we 
have, as in Art. (25), 

(period of moon about earth)* __ ma* . 
(period of sun about earth)* /wa s ^ ' 

therefore if a be properly chosen, 

2 _7n'a' 3 _ wV 8 
vn 3 z:4 4 
p,a ha 

47. If we examine the equations of Art. (31) it will be 
seen that we cannot solve the differentia 1 equation in u to 
the second order until we know the value of s accurately to 
that order; for, suppose 



where F^ contains the terms of the second order not yet 
found; and suppose that among the terms in F 2 there be 
on~ sin(Z0-X), then this value of $ substituted in the 
expression (D) gives 



= a (I - $i* sin 3 (#0 - 7) - 3&Z sin (g0 - 7) sin(?0-X)... 
and equation (A) becomes, so far as this term is concerned, 



[coB{fo + Z)0-( 7 + X)}].... 

Now, if it should happen that Ms a small quantity, then 
the two coefficients g -. / and g -f I will be both near unity, 
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and therefore these two terms, although of the third order 
in the differential equation, would become important in the 
value of u. Again, if I is near 2, then I - g will be near 
unity, and the same remark will apply. 

Such terms do not really occur in the value of $, but 
there is none the less a necessity for verifying the fact. 

If I is near unity, then the coefficient g I becomes small ; 
and although the corresponding term in u is still of the third 
order, it becomes of the second order, and therefore impor- 
tant in *, Art. (30). 

If I is not near 0, 1 or 2, the term will not affect, to 
that order, the values of u or of t. 

48. Again, it is necessary to have the value of u correct 

T 

to the second order, before we can compute r^-g to the fourth 

order, with a full assurance that we have secured all terms 
of that order in which the coefficient of the argument 
is nearly zero*; such terms becoming, as we have seen 
Art. (30), important in the value of t. 
For, assume that 



. T - . ,, _ , a 

then 77-. == ^TT-TT; - T\ - r 77rrr sm {(2 -2m) 0- 
AV 2A V {1 -f e cos (c0 - a) u J 



= - f m* (1...- 4 JF 2 ) sin {(2 - 2m) 6 - 2#j. 

Therefore if W 9 contain a term R cos (rd p) of the second 
order where r is near 2 (except 2 - 2wi) there will be a 



* This was pointed out by Mr. Walton in the Quarterly Journal of Mathematics, 
Vol. IX., p. 227; and in this edition the text had been modified in accordance 
with his suggestion j but the method given by him fo~ finding the value of u ia 
defective on account of its not mentioning the necessity for the previous deter- 
mination of *. 



44 LUNAR THEOKY. 

T 

corresponding term in ^- 3 of the fourth order with a small 

coefficient of 0, and therefore an important term in the 
value of t. 

49. For the sake of simplifying the expressions we shall 
generally omit a, $, and 7 in the arguments. If we remark 
that c and cc always appear together in the form c9 a, it 
will be sufficient to write cO l instead of c#-a; and, for 
the same reason, we may replace m6 + ft and gd y by 
m0 } and gO r The suffix (J will remind us of the omission, 
and we may at any stage re-introduce the omitted symbols. 



To compute s to the second order. 

50. Since we must first determine the value of s, it 
becomes necessary to examine the differential equation (B) 
to see whether any difficulty of a similar kind to those we 
have been here considering may arise, or whether the 
elements on which the computation depends are already 
known with sufficient accuracy. We fortunately find that 
all the terms of the right-hand member of the equation are 
of the third order, and therefore that the values already 
obtained will ensure accuracy to the second order, and we 
may safely begin our work here. 

51. To the first order we have 

s = k sin^, 

ds /IT /, (because q = 1 

-70 = A# COB00 = A;cosar0 t \ y 

d0 y * l y l (to first order, 

fa 
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- Ps 8 
In the expression for 77-;,- , which is of the third order, 

it will be sufficient to put a for %', and a for u^ then 



*lm*k sin (2 - 2m - 0) t , 

the other term of the third order being rejected because the 
coefficient of 6 is not near unity. 

T ds ds . 

The next term of equation (B\ is TO-S --,,;, and as -.--.is 
^ v J nu dv 7 dv 

T 
of the first order, it will be sufficient to have -, 3 correct 

to the second order, 

. im 1 Bin (2 -2111)^ 

fai - j e = - Sw** sin (2 - 2m) ^ COB^ 

= - %m?k sin (2 - 2m -^) 0^ 

neglecting the other term of the third order because the 
coefficient of 6 is not near unity. 

(7jj * * np 

-7L2 + S ) Ij*-*d9 

f T 
which =0 to the third order, since \j^ A d6 depending on 

the disturbing force is of the second order, and the other 
factor = to the first order. 



The equation in $ becomes 
d*s 



+ f m 9 A sin (2 - 2w - g) 0^ 



To solve it, assume 
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Substitute, and equate coefficient of like terma, 
k (1 - f) = - | l * ; .'. = (1 + }')* - 1 + fm f , 



^ {1 - (2 - 2 M - ,)<} = X* ; ... ^ = - 
Therefore, accurately, to the second order,* 



- 7) -f & sin {(2 - 2m -g}0- 2/3 + 7}... # 2 . 
The terra of the second order, having a coefficient near 
unity, need not be taken into account in finding the value 
of u (47) ; but it will have to be brought in subsequently 
to determine t. 

To determine u to the second order. 

52. We must compute the right-hand member of the 
equation (A), Art. (31), accurately to the second order, and 
include those terms of the third order which have a coefficient 
near unity. 

To the first order we have 



du . . . . (because c = 1 

-TTJ^ aecamcff t = -ae sin cv. 4 , 
dd l l (to first order, 

d'u 

d p + u = a. 

From the expressions (Z>) and (.E), Art. (31), we get 



* No complementary term P cos (0 - Q) is added ; for, though by the theory 
of differential equation? this would form a necessary part of the solution, we 
have seen Art. (39) that it cannot, in this shape, form a part of the correct value, 
hut will be comprised in the terms whose argument is g9 y. 



TO FIND t TO THE SECOND ORDER. 



47 



- f ) - 3e cosc0J {14-3 cos(2 - 2m) 0J, 
.1 ^ = 1 - f ' 4- ftf cos 2 9 e i - X - X cos (2 - 2m) 0, 

+ $m*e coseflj + -Jm'e cos (2 - 2m - c) X . 
The other terms of the third order, which arise in the 
development of the expression, are neglected because the 
coefficients of are not near unity, 
T 3m'a' 3 {! + <>' cos (wfl.-f)} 3 

_ _ __ _ I __ > __ 1 . 3 _iL a 

AV" aA "" 



= - | HI* [1 + 3e ; cos (771^ - f ) - 4c cosc^J sin (2 - 2m) l 
= - |m* sin (2 - 2m) l -f- 3?n s 6 sin (2 - 2m - c) lf 



i9 -;TI = |w 2 ae sin (2 2m) 0. si 
hu dv * ^ ti 

= f m*ae cos (2 - 2m - c) 0, ; 
rejecting the other terms of the third order according to (30). 



= | m 2 cos (2 - 2m) l - 3m 2 6 cos (2 - 2m - c) 0^ 

t^ [ 2 ,-^ 3 ^==3 7n 8 aC os(2~2m)0 l --6m 8 aecos(2--2m--c)0 1 . 
J j hu 



The equation in u becomes 



4- ym^e cos (2 - 2m - c) ^. 



To integrate this, assume 



u 
a 



1+A 
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Substitute, and equate the coefficients of like terms : 



C {I - (2 - 2w) 2 } = - 3w', 



Therefore 



# a) 
-] A; 2 cos 2 (^- 



U. 



V me cos {(2 - 2m - c) - 2/3 + a} 
This is the accurate value of u to the second, order 

To find t to the second order. 
53. The equation to be considered is 



and, making use of the values of s and u obtained to the 
second order, we must retrace our steps to recompute u and 

T 

j y - 8 so as to include those terms of the third order in u in 

/ i U> 

which the coefficient of the argument is near zero, and 

T 

similar terms to the fourth order in /2 - 3 Art. (48). To do 

ii> u 

this we need not repeat the whole work, but only so much 
of it as will enable us to pick out the additional terms in 
question. 
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-=1-3 {* sin .^ + %mk sin (2 - 2m - g) 0$ 

008(2-2/100,} 



in^ sin (2- 2m- 

cos (2 - 2m - 20) 1 - j]w<V cos (7*0, - f), 
r _ _ 3wV 3 (I + fi f COSW0, - )' si" (2 - 2w) t 
A V ~ 2/ 4 s a 4 ( L + e cos c0 t - -}/,;' cos 2y ^ ) 4 " 

= - f m* sin (2 - 2?w) l (I + e cos^^ - {fc* cos2^^} " 4 
= - ?/i s sin (2 - 2i) 0, {(I + c COSC0,)' 1 4- &* COS2//0J 
= ~ -Jjn* sin (2 -2/w) 0, {...10e 
= - ?>A 3 sin (2 - 2)/0 0, {...o^ 
= Y J wVsin(2-2wt-26)0 l ~2^^ sin(2-2?-2^) r 
These are the only important additional terms. 

r T 

1 72 3 ^0=- V we* cos (2- 2m -2c) 0,- ^tnlf co(2-2fi-2y) lf 

because 2 -2w - 2c = - 2m + } w 2 , and 2 - 2?/z - 2g = - 2m f m 2 . 
Substituting in the differential equation for w, the new 

terms in y^ -77. will be of the fourth order and must bo 
li u au 

neglected, but 2 (~i^ + u ) lit *<$ wil1 bring in the two 
terms 

- *me*a cos (2 - 2m - 2c) l and - %mk*a cos (2 - 2m - 2/7) l ; 
and we get, so far as important additional terms are concerned, 

a 

+ (- A + 1) w ^ cos ( 2 - 2m - 2 ,7) e r~ iwY cos 
and the corresponding terms in u will also be of the third 
order, viz. 

- =. . . ^ me* cos (2 - 2m - 2c) 0, -f ^i* a cos (2 - 2m - 2g) 0, 

- |mV cos(?n0 t - f). 
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The complete value of u therefore, so far as important terms 
are concerned for determining the value of , is 



- = 1 _ | A 3 - \m* + e cosc0 l - 



+ m a cos(2 - 2m) l + igme cos(2 - 2m - c) l 

H- *' c <> s (2 - 2m - 2c) l + -&iwS" cos (2 - 2m - 20) X 



mY 



and the complete value of ly^s ^ W ^ 

j Afc ZA 

/T 
7T~8 ^= f 8 cos (2 - 2m) ^ - V^* cos (2 - 2m - 2c) t 
ft, U 

- tyntf cos (2 - 2m - 2#) ^, 

the other term of the third order has no influence on tj the 
coefficient of not being small. 

54. We are now in a position to find the value of t to 
the second order : 



If - = 1 -f e cosc0 t -f W where W contains all terms of the 
a 1 

second and third orders in the expression above, then 



, {1 -2 (e COSC0.+ JF) + 3 ( e 
-^ 2 (1 - 2e COSC0, - 2 W+ 3e 8 cos a 
^ (i 4. | e s - 2^ cosc^- 2 IT+ f e 2 
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The only part of 60 TFcoscfl, which gives a term of the third 
order in which the coefficient of the argument is small, is 

Gecoscflj {...-h *gme cos(2- 2m c)^ ...}, 
of which one term is 

cos (2 - 2m - 2c) r 



Therefore ^-5 is equal to 



!1 + 1* 1 -f f If + m'-2e cosc^ + f e* 
+ } 2 k* cos2g0 t - 2m 2 cos (2 - 2m) t 
- ^me cos (2 - 2m - c) 0, -f 3m V cos(m^ - ) 
- ^me* cos (2 - 2m - 2c) ^ fwtfi* cos (2 20i - 
Also from above 



fT 
JhV 



cos (2 - 2m - 2g) 0,. 

We have now to multiply these results together, and we see 
that the terms having for coefficients *gme* and fm/c" will 
disappear in the product.* These were the terms which 

T 
were originally of the fourth order in 7.^-3 . 

li U 



(1 4 
4 
- 



4 * coB2jr^ - V w C03 ( 2 - 

cos (2 - 2m - c) 0, + 3mV cos(rw^ - ?) 



* The principle of the superposition of small disturbances might have led us 
& priori to anticipate the disappearance of these terms: -for, when approximating 
to the second order we may, according to this principle, calculate separately the 
effect of the disturbing action of the sun, which is of the second order, and that 
of any other independent disturbance of the same or of a higher order, and then 
add the results. No term therefore which bears a trace of the action of both 
these forces should present itself to the second order. 

Now since neither k nor e enters into the argument of any term, it is obvious 
that a term whose coefficient contains A; 2 or e 2 at any stage will retain them 
through all subsequent operations ; and if the terms having (2 2wi 2^) l 
or (2 - 2/ - 20) 0! for argument could remain in the result to the second order 
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Let li + e " + 



sin {(2- 2m) -2/3} 



therefore p = Tic? (1 \ f 7c 2 m 2 ) to the third order ; 
therefore, multiplying by p and integrating, we get, still to 
the second order, 



.... 



no constant is added, the time being reckoned from the 
instant when the mean value of 6 vanishes, for the reasons 
explained in (Art. 36). 

Coordinates expressed in terms of the time. 

55. The preceding equations f , $,, 2 give the reci- 
procal of the radius vector, the latitude and the time in 
terms of the true longitude ; but the principal object of the 
analytical investigations of the Lunar Theory being the 
formation of tables which give the coordinates of the moon 
at stated times, we must express w, s, and 6 in terms of t. 

To do this, we must reverse the series pt = ff &c., and 
then substitute the value of in the expressions for u and s. 

Now =pt H- 2e sin (c9 - a) to the first order 

sm(cpt a) ................... ; 



they would there liavc IP or ^ respectively in the coefficient, since they have 
them when they first appear. Bnt X, 2 , e 2 and m being independent, such terms 
must have arisen from the combination of others which originally (in the 
expressions of the forces) involved these quantities separately : terms like 
k z eos'2(#t) y) and e- cos 2 (rO a) combined with terms involving mi that 
is, terms which are representatives of disturbing forces of the second order (see 
note p. 31) combined wi+h others depending on the sun's action. Such terms 
wust therefore be of an order beyond the second. 
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therefore cO a = cpt a + 2e sin (cpt a) to the first order, 
2e sin (c0 a) = 2e {sin (cp a) + 2e sin (cpt a) cos (cp a)} 

to the second order? 
= 2esir(c^-a) + 2e 2 sin2(c/^- a) ............ ; 

and as and pt differ by a quantity of the first order, they 
may be used indiscrirnirately in terms of the second order; 
therefore 

0=jpt+2e sin (cpt ~ a)-f- -|e* sin 2 (cpt a) 
-Jr/v 2 sin 2(^-7) 
+ V m* sin {(2 - 2m) pt - 2/3} 



3?e' sin (nipt + /8 - 



... 0' a . 



56. In the value of u given in Art. (52), substitute pt for 
in terms of the second order, and pt + 2e sin(cpt - a) in the 
term of the first order ; then 



^A^ 3 cos 2 (gpt - 7} 
+ m 2 cos {(2 - 



57. Similarly, the expression for s becomes 
s = A sin {(<7p 7) 4 2<s sin (cp - a)} 

+ fynk sin {(2 - 2m - <j) pt - 2/3 + 7} 



or 



s = 



sn t- 



- e sin {(</ c)^)^ + a 7} 
. 4- m sin {(2 - 2m - 0) /tf - 2/3 + 7} 
The expression for 5 is more complex in th ; s form than when 
given in terms of the true longitude 0. 
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MoorCs Parallax* 

58. If P be the moon's mean parallax, and FT the parallax 
at the time t, 

_ radius of earth R - , . ~ ,. , , 

n= __ -- - ~ = = jBn(i -is 8 ) to the third order, 

distance of J 1 .. x 

e= JBw {1 J& 2 -f J/c 2 cos2 (gpt 7)} to the second order, 

11 - k 3 \m* e 2 -I- e cos (cpf a) + e* cos2 (qp - a) 
f w 2 cos {(2 - 2m) pi - 2/3} 
H- ^ 7?ze cos {(2 - 2m - c}pt - 2/3 -f a} ; 
but jP= the portion which is independent of periodical terms, 



1 -f e cos(cpt-a) -f e 8 cos 2 (cpt a) 
therefore IT = P , + m* cos {(2 - 2?w) pt - 2/8} 

-f ^ 5 wie cos {(2 - 2m - c) ^pf - 2/3 -f a} 

neglecting terms of the third order. We see that, to the 
second order, the variable part of the parallax is independent 
of the inclination. 

59. Here we terminate our approximations to the values 
of w, , and 0. If we wished to carry them to the third 
order, it would be necessary to include some terms of the 
fourth and fifth orders according to Art .(30), and the ap- 
proximate values of P, T, and $, given in Art. (24), would 
no longer be sufficiently accurate, but we should have to 
recur to the exact values, and from them obtain terms of 
an order beyond those already employed.* The process 
followed in the preceding pages is a sufficient clue to what 
would have to be done for a higher approximation. 

* See Paralfactic Inequality, Art. (105). 
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The coordinates u' and & of the sun's position are, by 
the theory of elliptic motion, known in terms of the time , 
and t is given in terms of the longitude by the equation 
G 2 . Hence u' and & can be obtained in terms of 0; but 
it will be necessary to take into account the slow progressive 
motion of the sun's perigee, which we have hitherto neglected* 
This we may do by writing cff for ff f, c' being & 
quantity which differs very little from unity.* 

These values of u 1 ^ 6 ', together with those of u and $ 
given by U^ and 8^ are then to be substituted in the 
corrected values of the forces, and thence in the differential 
equations. The integrations being performed as before will 
give the values of u, s, and t in terms of to the third 
order, and from these, as in Arts. (55), (56), and (57), may bo- 
obtained w, a, and 6 in terms of t. 

60. More approximate values of c and g are obtained at 
the same time, by means of the coefficients of cos (c0 a} 
and sin (gd 7) in the differential equations, (see Appendix,, 
Arts. 102 and 103). 

61. The values to the fourth order are then obtained 
from those to the third by continuing the same process, and 

* * En reflechissant sur les termes que doivent introduire toutcs Ics quantitca 

* pr6cedentes, on voit qu'il sc peut glisser des cosinus de Tangle dont nous avons 
*vu le dangereux effet d'amener dans la valeur de u des arcs au 1'su de leura 

* cosinus j de tcls termes viendront, par exemple, de la combinaison des cosinus 
'de (1 m) avec des cosinus de mQ 

* Pour eviter cet inconvenient qui oterait 3. la solution precedento 

'Tavantage de convenir & un aussi grand nombre de revolutions qu''on voudrait, 
4 et la priverait de la simplicite et de 1' universalite si precieuses en matliematiques, 
*il faut commencer par en chercher la cause. Or, on decouvre facilement que 
'ces termes ne viennent que de ce qu'on a suppose fixe 1'apogee du soleil, ce 
' qui n'est pas permis en toute rigueur, puisque quelque petite que soit sur cet astre 

* 1'action de la lune, elle n'en est pas moins reelle et doit loi produire un mouvement 
' d'apogee quoique trds lent la verit'. Clairaut, Theorfc de la Lune, p, 55. 
2me Edition. 
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so on to the fifth and higher orders ; but the calculations are 
so complex that the approximations have not been carried 
beyond the fifth order, and already the value of 6 in terms 
of t contains 128 periodical terms, without including those 
due to the disturbances produced by tlie planets. The coeffi- 
cients of these periodical terms are functions of ?n, 0, e', , 

c, #, 7c, and are themselves very complicated under their 
literal forms : that of the term whose argument is twice the 
difference of the longitude of the sun and moon, for instance, 
is itself composed of 46 terms, combinations of the preceding 
constants. 

See Pontdcoulant, Systime du Monde^ torn. IV., p. 572. 



CHAPTER V. 

NUMERICAL VALUES OF THE COEFFICIENTS. 

62. Having thus, from theory, obtained the form of the 
developments of the coordinates of the moon's position at 
any time, the next necessary step is the determination of the 
numerical values of the coefficients of the several terms. 

We here give three different methods which may be 
employed for that purpose, and these may, moreover, be 
combined according to circumstances. 

63. First Method. The values of the constants p, m, a, 
PI 7, f which enter into the arguments, and of the additional 
ones a, 7c, 0, &c. which enter into the coefficients^ of the 
terms in the previous developments, may be obtained with 
great accuracy from observation. . 

For this we must employ observations separated by very 
long intervals, such, for instance, as ancient and modern 
eclipses; also particular observations of the sun and moon 
made when the bodies occupy certain selected positions; 
and other observations of a special character. 

When the values of the elements have been so obtained, 
then the theoretical values of the coefficients may be com- 
puted by substitution in the analytical expressions. 

64. Second Method. Let the constants which enter into 
the arguments be determined as in the first method ; and let 
a large number of observations be made, from each of which 
a value of the true longitude, latitude, or parallax is obtaine' 1 , 
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together with the corresponding value of t reckoned from 
the fixed epoch when the mean longitude is zero. Let these 
corresponding values be substituted in the equations, each 
observation thus giving rise to a relation between the unknown 
constant coefficients. 

A very great number of equations having been thus 
obtained, let them, by the method of least squares or some 
analogous process, be reduced to as many as there are 
coefficients to be determined. The solution of these simple 
equations will give the required values. 

This method, however, would scarcely be practicable in 
a high order of approximation. For instance, in the fifth 
order, as stated in Art. (61), each observation would give 
rise to an equation containing 129 unknown quantities and 
the immense number of equations so obtained would have 
to be reduced to 129 equations of 130 terms each. 

65. Third Method. When the constants which enter into 
the arguments have been determined by the first method, we 
may obtain any one of the coefficients independently of all 
tlw others by the following process, provided the number of 
observations be very great. 

Let the form of the function be 



and let it be required to determine the constants -4, J5, 0, &c. 
separately; 0, <, &c. 3 being known functions of the time. 

Let the results of a great number of observations corre- 
sponding to values t , 8 , g , &c., < 17 2 , < s , &c., be F 1} F a , F 3f 
&c. ; so that 

FjS A -f-Bsin^-f Csin^-f &c., 
F 2 = A + B sin 2 4- tfsin^ + &c., 



V '=-A -f B sin^-h <7sin<f> 4 &c. 
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Now, n being very great, we may assume that the sum of 
the positive values of each periodical term will be about 
counterbalanced by the sum of its negative values; and 
therefore, that if we add all the equations together these 
terms will disappear ; 

., , , 

therefore A 



which determines the non-periodic part of the function, 

To determine B. Let the observations be divided into 
two sets separating the positive and negative values of sin ; 
then the other periodical terms, not having the same period, 
may be considered as cancelling themselves in adding up the 
terms of each set. Let there be r terms in the first set 

and s terras in the second, and let F', F", F r be the 

values of F corresponding to positive values of sin 0, which 
values we may assume to be uniformly distributed from 

sin0 to sin-Tr, and therefore to be sin 80, sin 280, sin r. 80, 

where r.80 = 7r. 

And, again, let F,, F /y , F //y , F,, be the values of F 

corresponding to the negative values of sin0 ; viz., sin (- A0), 

sin(-2A0), sin(-s.A0), where s.A0 = 7r. Then, 

F' = -4 + J9sinS0 -f Csin</>' +..., F y = A B sin A0 + (7sin^ y 4..., 



F r = 
therefore 






therefore 



B [ T 

**]* 



7T 



therefore 



therefore 



A 
J ; 
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o, r T? * ( v ' + V"+...+ V' F,+ F,, '+...+ F 
therefore 5 = - ^ -- - 

and in a similar manner may each of the coefficients be 
independently determined.* 

66. When the periods of two of the terms differ bu* 
slightly for instance if Q and < go through their periodic 
variations very nearly in the same time, the method could 
not then with safety be applied ; for, since the same values 
of 6 and would very nearly recur together during a longer 
time than .that through which the observations would ex- 
tend, the two terms would be so blended in the value of 
V that they would enter nearly as one term the difference 
between and <j> would be very nearly the same at the 
end as at the beginning of the series of observations. 

* If r and s are not sufficiently great to allow us to substitute f 9 v Bin Odd for 
E^ siu Q . <50, we must proceed as follows : 

V'+ V" 




_. .. . r, + V H + ...... + r s . 2B / 7r 2 \ 

Similarly, -* - "-^ - = A - ^1 - T V^J J 

TT ,y+V"+ ...... + pr ^ y f + r y/+ ...... + yx 

^ 1 r ^ r 5 ' 



therefore 

~ 
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67. Let us suppose the periods to be actually identical, 
so that < = -f a, a being some constant angle ; then 

B sin 6 -f C sin <j> 

may be written (S+C ^os a) sin 6 -f C sin a cos 0, 
T V=A + (J?+(7cosa) sin0+(7sinacos0-f ....... 

If now we divide tLe observations, as before, into two 
sets, corresponding to the positive and negative values of 
sin#, the terms involving cos# will disappear in the sum- 
mation of each set ; and, following the process of the method, 
we shall find the value of 

B -f C cos a = M suppose. 

Dividing again into two sets corresponding to the positive 
and negative values of cos#, the terms in sin0 t will be can- 
celled, and the same process will give the value of 

C sin a = N suppose. 

Treating the observations in the same way with respect 
to the angle $, we get two results, 



B 
from these four equations we easily get 



M'N+MN' 

*-" '- 7 * 7 



'j N' are connected by the equation of condition, 



68. When the periods of and <j> are nearly, but not 
exactly, the same, this equation of condition will not hold, 
and the preceding values of B and C would not be exactly 
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correct ; but yet they would be very approximate, especially 
if the mean between the two values of B be taken. 

Or, ' we may, after having taken one of these slightly 
erroneous values for J?, make a further correction by estab- 
lishing as it were a counterbalancing error in the value of 
C. Let 2? be the value so found for B; then, from the 
F of each of the observations subtract the value B' sin 0, 
the result ?7wiU be very nearly equal to A + C sin ^ + &c., 
and from the n equations 

U^A + 



a value (7 of C will be obtained, by the rule, which will be 
very approximate, and, at the same time, agree better with 
J3' in satisfying the equations than C itself would do. 

69. When two terms whose periods are nearly equal do 
occur, it is plain, by examining the values of M and M ', that 
the errors which would be committed by following the rule, 
without taking account of this peculiarity, would be the taking 
JB+ G cos a and(7-f 13 cos a for B and C respectively. 



CHAPTER VI. 



PHYSICAL INTERPRETATION. 

70. The solution of the problem which is the object of 
the Lunar Theory may now be considered as effected ; that 
is, we have obtained equations which enable us to assign 
the moon's position in the heavens at any given timo to 
the second order of approximation ; we have explained how 
the numerical values of the coefficients in these equations 
may be determined from observation ; and we have, more- 
over, shewn how to proceed in order to obtain a higher 
approximation.* 

It will, however, be interesting to discuss the results we 
have arrived at : to see whether they will enable us to form 
some idea of the nature of the moon's complex motion ; also 
how far they will explain those inequalities or departures 
from uniform circular motion which ancient astronomers had 
observed ; but which, until the time of Newton, were so 
many unconnected phenomena; or, at least, had only such 
arbitrary connexions as the astronomers chose to asrign, by 
grafting one eccentric or epicycle on another as each newly 
discovered inequality seemed to render it necessary. 

It is true that our expressions, composed of periodic 
terms, are nothing more than translations into analytical 

* The means of taking into account the ellipsoidal figure of the earth and the 
disturbances produced by the planets, are too complex t ) form part of an intro- 
ductory treatise. For information on these points reference may be made to- 
Airy'a Figure of the Earth. Pontecoulant's Systeme du Monde, vol. iv. 
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language of the epicycles of the ancient;* but they arc 
evolved directly from the fundamental laws of force and 
motion, and as many new terms as we please may be ob- 
tained by carrying on the same process ; whereas the epicycles 
of Ilipparchus and his followers were the result of numerous 
and laborious observations and comparisons of observations , 
each epicycle being introduced to correct its predecessor 
when this one was found inadequate to give the position 
of the body at all times: just as with us, the terms of the 
second order correct the rough results given by those of 
the first ; the terms of the third order correct those of the 
second, and so on. But it is impossible to conceive that 
observation alone could have detected all those minute 
irregularities which theory makes known to us in the terms 
of the third and higher orders, even supposing our instruments 
far more perfect than they are ; and it will always be a 
subject of admiration and surprise, that Tycho, Kepler, and 
their predecessors should have been able to feel their way 
so far among the Lunar inequalities, with the means of 
observation they possessed. 

LONGITUDE OF THE MOON. 

71. We shall first discuss the expression for the moon's 
longitude, as found Art. (55), 

6 -=pt 4- 2c sin (cpt - a) -f ^ sin 2 (cpt - a) 
+ Y' me sin {("2 - 2m - c) pt - 2/3 -f a} 
+ V f sin {(2 - 2m) pt - 2/3J 



The mean value of is pt] and in order to judge of the 
effect of any of the small terms, we may consider them 



* See Whewcli's 1/htnry <>/* tht /</<///< 
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one at a time as a correction on this mean value pt, or 
we may select a combination of two or more to form this 
correction. 

We shall have instances of combination in the elliptic 
inequality and the evcUion 9 Arts. (73) and (77); but in the 
remaining inequalities each term of the expression will form 
a correction to be considered by itself. 

72. Neglecting all the periodical terms, we have 

6-fr 
(10 
lit = *> 

which indicates uniform angular velocity; and as, to the 
same order, the value of u is constant, the two together 
indicate that the moon moves uniformly in a circle, the 

2?r 
period of a revolution being - , which is, therefore, the 

expression for a mean sidereal month, or about 27 days.* 
The value of p is, according to Art. (54), given by 



and as m is due to the disturbing action of the sun, we see 
that the mean angular velocity is less, and therefore the mean 
periodic time greater than if there were no disturbance. 

Elliptic Inequality or Equation of the Centre. 

73. We shall next consider the effect of the first three 
terras together : the effect of the second alone, as a correction 
tj will be discussed in the Historical Chapter, Art. (115). 



* The accurate value was 27d. 7h. 43m. ll*2tiis. in the year 1801. See 
Art. (10*). 

K 
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6 pi -f 2e sin (cpt - a) 4- \e* sin 2 (cpt - a), 
which may be written 
6=pt+%e sin [pt~ {a. -f-(l - c)pt}]+&e 2 sin 2[p$ - {+ (1 

But the connexion between the longtitude and the time in 
an ellipse described about a centre of force in the focus i fl / 
Art. (13), to the second order of small quantities: 

= nt + 2e sin (nt - a') + \<? sin 2 (nt - a'), 

where n is the mean motion, e the eccentricity, and a' the 
longitude of the apse.* 

Ilence, the terms we are now considering indicate motion 
in an ellipse; the mean motion being p, the eccentricity e, 
and the longitude of the apse a 4(1 c) pt; that is, the 
apse is not stationary but has a progressive motion in longi- 
tude, uniform, and equal to (1 - c)p. 

74. The two terms 2e sin (cpt - a) + f e* sin 2 (cpt - a) con- 
stitute the elliptic inequality , and their effect may be further 
illustrated by means of a diagram. 

Let the full line AMB re- 
present the moon's orbit about /^'"^ ^ X 
the earth E, when the time t ~^ ^ * "" " 5 
commences, that is, when the 
moon's mean place is in the prime 
radius JET, from which the lon- 
gitudes are reckoned. 

The angle Y2L4, the longi- 
tude of the apse, is then a. At the time , when the moon's 
mean longitude is VEM=pt, the apse line will have moved 



* The epoch which appears In the expression of Art. (13) is here omitted j 
a proper assumption for the origin of t, as explained in Art. (36), enabling us 
to avoid the &. 
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in the same direction through the angle AEA' = (1 c) TJEJf, 
and the orbit will have taken the position indicated by the 
dotted ellipse. The true place of the moon in this orbit, 
so far as these two terms are concerned, will be ?H, where 

MEm = 26 sin (cpt - a) + f e* sin 2 (cpt - a) 
2e sit A' EM 4 e' sin 2 A 1 EM 



which, since e is about ^, is positive from perigee to apogee, 
and therefore the true place before the mean; and tho 
contrary from apogee to perigee : at the apses the places 
will coincide. 

75. The angular velocity of the apse is (1 c) p, or, if 
for c we put the value found in Art. (52), the velocity will 
be $m*p. Hence, while the moon describes 360", the apse 
should describe f wi*.3GO= 1 nearly, mi being about ^. 

But Hipparchus had found, and all modern observations 
confirm his result, that the motion of the apse is about 3 in 
each revolution of the moon. See Art. (118). 

This difference arises from our value of c not being repre- 
sented with sufficient accuracy by 1 w 2 . 

Newton himself was aware of this apparent discrepancy 
between his theory and observation ; and we are led, by his 
own expressions (Scholium to Prop. 85, lib. in. in the first 
edition of the Principia), to conclude that he had got over 
the difficulty. This is rendered highly probable when we 
consider that he had solved a somewhat similar problem 
in the case of the node ; but he has nowhere given a state- 
ment of his method : and Clairaut, to whom we are indebted 
for the solution, was on the point of publishing a new hy- 
pothesis of the laws of attraction, in order to account for 
it, when it occurred to him to carry the approximations +o 
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the third order, and he there found the next term in the 
value of c nearly as considerable as the one already obtained. 
See Appendix. The value is 



.-. I - c = lm* -f V 
/. (1 - c) 360 = (1 4- iW) ( v alc round previously) 
= 22" nearly, 

thus reconciling theory and observation, and removing what 
Lad proved a great stumbling-block in the way of all as- 
tronomers.** 

When the value of c is carried to higher orders of ap- 
proximation, the most perfect agreement is obtained. 

The motion of the apse line is coiibidcred by Newton 
in his Principle!,) lib. I., Prop. GG, Cor. 7. 

Emotion. 

70. The next term + \^m sin {(2 - "2>n - o) pi - 2/3 + a} 
in the value of 6 has been named the Evectlon. We shall 
consider its effect in two different ways. 

First, by itself, as forming a correction on pt^ 

=zpt + me sin {(2 - 2m - c] pt - 2/3 + a}. 
Let J = pt = moon's mean longitude at time ^ 
^ =5 mpt + & = sun's ................................. j 

a r = (1 - c) pt 4 a = mean longitude of apse ......... , 

then 



V s w ^ sin t 2 {^ - 

V;;w sin [2 (j - ) _ (5 - a 1 )}. 



See Dr. WlicweU's llri'li/wahr Tnatiw, 
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The effect of this term will therefore be as follows : 
In syzygies 

O^pt ^mc sin (5- a'); 

or the true place of the moon will be before or behind the 
mean, according as the moon, at the same time, is between 
apogee and perigee or between perigee and apogee. 
In quadratures 

^pt -f- if me sin (J) a'), 

and the circumstances will be exactly reversed. 

In both cases, the correction will vanish when the apse 
happens to be in syzygy or quadrature with regard to the 
sun at the same time as the moon. 

In intermediate positions, the nature of the correction 
is more complex, but it will always vanish when the sun 
is at the middle point between the moon and the apse, 
or when distant 90 or 180 from that point; for if 

G = ~^"~ - r . 90, where r = 0, 1, or 2, 

sin [2 (J - o) - (3) - a')] = sin () + a' - 20) 
= sin r. 180 
= 0. 

77. The other and more usual method of considering 
the effect of this term is in combination with the two terms 
of the elliptic inequality, as follows : 

To determine the change in the position of the apse and in 
the eccentricity of the moorfs orbit produced ly the evectwn. 

Taking the elliptic inequality and the evection together, 
we have 

0~pt + 2e sin (cpt - a) 4- 1<? sin 2 (cpt - i) 

+ * me sin {(2 - 2m - c]pt - 2/3 -f a}. 
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Let a' be the longitude of the apse at time t on supposition 
of uniform progression, 
O ................ .............. sun ............ ; 

whence a! = (I c}pt + a, 

O = mpt + 0. 

And the ahove may be written t 

6 =*pt 4- 2e sin (cpt a) -f \<? sin 2 (r/; a) 

-f- l l>mc sin {qrf - a 4- 2 (a' - 0)} ; 
and the second and fourth terms may be combined into one, 



if we assume E cos 8 = e + J b 5 me cos 2 (a' o )> 
and E sin 8 = *gme sin 2 (a' ) ; 

whence tan 8 and E may be found ; and approximately, 
E=e{\ + 1 t fmcos2 (a 1 o)}, 
8 = ^w* sin 2 (a' o). 

The term Jc* sin2(<2rt-a) may also, therefore, to the 
second order, be expressed by 



and the longitude becomes 

6 ~pt + 2E sin (ept - a -I- 8) + f E* sin 2 (cp* - a + 8), 
or = JB -f 2^sin (^ - a + 8) + $J?* sin 2 ( jtf - a' + 8) ; 
but the last two terms constitute elliptic inequality in an orbit 
whose eccentricity is E and longitude of the apse a' 8 ; 
therefore the evcction, taken in conjunction with elliptic 
inequality, has the effect of rendering the eccentricity of the 
moon's orbit variable, increasing it by l ^ me when the apse- 
line is in syzygy, and diminishing it by the same quantity 
when the apse-line is in quadrature ; the general expression 
for the increment being *gme cos 2 (a' o). 
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And another effect of this terra is, to diminish the longitude 
of the apse, calculated on the supposition of its uniform 
progression, by the quantity 8 = y> m sin 2 (a' - 0) ; so that 
the apse is behind its mean place when in the first and third 
quadrant in advance of the sun, and before its mean place 
in the second and fourth. 

The cycle of these jhangcs will evidently be completed 
in the period of half a revolution of the sun with respect 
to the apse, or in about -/<. of a year.* 

78. The period of the evcction itself, considered inde- 
pendently of its effect on the orbit, is the time in which 
the argument (2 2m c)jjt 2/3 + a will increase by 2?r. 

Therefore period of evection 

__ 2-7T __ mean sidereal month 

(2 2m c] p 2 2/>i c 

mean sidereal month 27 \ days 
-- --- 1 



= 31^ days, nearly, f 

Newton has considered the evection, so far as it arhca 
from the central disturbing force, in Prop. G6, Cor. 9, of 
the Principia. 

Variation. 

79. To explain the physical meaning of the term 
sin {(2 - 2wi)^ - 2/3}, 



* The change of eccentricity and the variation in the motion of the 
follow the same law as the abbcfcsa and ordinate of an ellipse referred to 
its centre : for if e = x and *5 #, then 



t The accurate value is ol '8119 days. 
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in the expression for the moon's longitude, 

=pt + y in 2 sin {(2 - 2m) pt - 2/3}. 
Let }) represent the moon's mean longitude at time t 9 
.................. sun's 

therefore )) =pt, 



and the value of 6 becomes 

0=^4- y 3 sin2(j- 0), 

which shews that from syzygy to quadrature the moon's 
true place is before the mean, and behind it from quad- 
rature to syzygy; the maximum difference being Vwi 2 in 
the octants. 

The angular velocity of the moon, so far as this term 

is concerned, is 
7/a 

i/ (1 - m) vi* p cos2 (5 - O), 



-P i 1 + V 111 * cos2 (5 ~" )h nearly, 

which exceeds the mean angular velocity p at syzygics, is 
equal to it in the octants, and less in the quadratures. 

This inequality has been called the Variation, its period 
is the time in which the argument (2 2m) pt 2^8 will 
increase by 2ir ; 

. , ~ . .. 2-7T mean synodical month 

/. period of variation = -~r = - - - - 
1 (2 - 2m) p 2 

= 142 days, nearly.* 

80. The quantity y m 9 is only the first term of an endless 
scries which constitutes the coefficient of the variation, the 
other terms being obtained by carrying the approximation 

* The accurate value is U* 765*294 days. 
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to a higher order. It is then found that the next term in 
the coefficient is J-j!/M 3 , which is about -fa of the first term ; 
and as there are several other important terms, it is only 
by carrying the approximation to a higher order (the 5th 
at least) that the value of this coefficient can be obtained 
with sufficient accuracy^ from theory. In fact, VM* would 
give a coefficient of 2G' 27" only; whereas the accurate value 
is found to be 39' 80". 

The same remark applies also to the coefficients of all 
the other terms. 

81. As far as terms of the second order, the coefficient 
of the variation is independent of e the eccentricity, and k 
the inclination of the orbit. It would therefore be the samo 
in an orbit originally circular, whose plane coincided with 
the plane of the ecliptic: it is thus that Newton has con- 
sidered it. Princip. prop. GO, cor. 3, 4, and f>. 

Annual liquation. 

82. To explain the jdty steal meaning of the term 

3w' sin (nipt -f - ) 
in the expression for the moon's lonyitwle. 
0=2>t-- Sme' sin (mpt -f - f), 

=pt3mc' sin (longitude of sun - longitude of sun's perigee), 
= pt - 3me sin (sun's anomaly). 

Hence, while the sun moves from his perigee to his 
apogee, the true place of the moon will be behind the mean ; 
and from apogee to perigee, before it. The period being an 
anomalistic year, the effect is called Annual Equation. 
Differentiating 6 we get 

rlQ 

T ~P ft "~ 3//tV cos (sun's anomaly)}. 

eft 
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Hence, so far as this inequality is concerned, the moon's 
angular velocity is least when the sun is in perigee," that 
is at present about the 1st of January, and greatest when 
the sun is in apogee, or about the 1st of July. 

The annual equation is, to this order, independent of the 
eccentricity and inclination of the moon's orbit, and therefore, 
like the variation, would be the same in an orbit originally 
circular. Vide Newton, Primipia, prop. 66, cor. 6. 

Reduction. 

83. Before considering the effect of the term 

A' 51 
-~ 

which, as wo shall see Art. (89), is very nearly equal to 
the difference between the longitude in the orbit and the 
longitude in the ecliptic, it will be convenient to examine 
the expression for the latitude of the moon, and to see how 
the motion of the node is connected with the value of g. 

LATITUDE OF THE MOON. 

84. The expression found for the tangent of the latitude,* 
Art. (51), is 

s = Jc sin (gO - 7) 4 mk sin { (2 - 2m - g) - 2/8 + 7}. 
If we reject all small terms, we have 
8 = 0, 

or the orbit of the moon coinciding with the ecliptic, which 
is a first rough approximation to its true position. 



* This expression for the tangent of the latitude is more convenient than 
that which gives it in terms of the mean longitude, Art. (57), on account of 
the "ess number of terms involved. See Pontecoular% vol. iv. p. 030. 
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85. Taking the first term of the expansion 

$ = isin(j/0 7), 
we may write it 




Let TNm be the ecliptic, N the 
moon's node when her true longitude 
is 2ero, and let M be the position 
of the moon at time , m her place 
referred to the ecliptic ; 

therefore T-W*7, Tw = 0, tanJI///t = s. 

Take NN'=z(g l}0 in a retrograde direction; and join 
MN f by an arc of great circle ; 

then tan.i/#i= tanl/JV / wisinJV r 't, 

or s = tan MN'm sin [0 - {7 - (g - 1) 0}] ; 

which, compared with the value of s given above, shews that 
MN'm = tan" 1 /; is constant, and therefore the term k sin (y0-~y) 
indicates that the moon moves in an orbit inclined at an 
angle tan' 1 ^ to the ecliptic, and whose node regredes along 

the ecliptic with the velocity (# l)-y ) or with a mfian 
velocity (g - 1) p. 

86. Hence the period of a revolution of the nodes 
2?r one 8l d erea l month 



but, from Art. (51), the value of g = 1 + w* ; 

, - i /? i A - * i one sidereal month 

therefore period of revolution of nodes = -- ,._ - 

= 6511 days, nearly. 

This will, for the same reason as in the case of the apse, 
Art. (75), be modified when we carry the approximation 
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to a higher degree ; this value of g is, however, much more 
accurate than the corresponding value of c, for the third 
term of g is small ; the value to the third order being (see 
Appendix, Art. 102) 

, , i r i 4- c 1 i one sidereal month 
and the period or revolution ot the nores= -~- :} - 3 a 

= 6705 days, nearly. 

This is not far from the accurate value as given by 
observation, and when the approximation to the value of (j 
is carried to a higher order, the agreement is nearly perfect. 

The true value is 6793'39 days, that is about ISyrs. 7 nio. 

Ejection in Latitude. 

87. To explain the variation of the inclination and the 
irregularity in the motion of the node expressed % the term 
+ link sin {(2 - 2m - </) - 2/3 -f 7}. 

This term, as a correction on the preceding, is analogous 
to the evection as a correction on the elliptic inequality. 
Taking the two terms together, 

s = k sin ( g6 - 7) + g w* sin {(2 - 2m -<J)6- 2/8 + 7}. 
Let ]) = longitude of moon = 0, 

O= sun =m0 + }, 

2= node =7~(^ 1)0; 

therefore s = k sin (5 - ft ) i- *mk sin {5 - ft - 2 (O - ; 
Now these two terms may be combined into one, 

a = A" sin (5- Q -A), 
if K cos A = k + $mk cos 2 (O Q), 

A' sin A = *mk sin 2 ( ^ Q ), 
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whence A and K may be found ; and, approximately, 



-S3), 
but the equation 

/? = JTsinfJ V Q A) 

represents motion in an orbit inclined at an angle tan^JT to. 
the ecliptic, and the longitude of whose node is Q -f A. 

This term has therefore the following effects : 

1st. The inclination of the moon's orbit is variable, it 
tangent increases by -*mk when the nodes are in syssygies, 
relatively to the sun, and decreases by the same quantity 
when they are in quadrature ; the general expression for 
the increase being *mk cos2(0 Q). 

2nd. The longitude of the node, calculated on supposition 
of a uniform regression, is increased by A = jjw sin 2 ( w ), 
so that the node is before its mean place while in the first 
or third quadrant in front of the sun, and behind it in tho 
second and fourth. Priiicipia^ book lir., props. 33 and 35. 

The cycle of these changes will be completed in tho 
period of half a revolution of the sun with respect to the 
node, that is, in 173*21 days, not quite half-a-year. 

88. The tangent of the latitude has here been obtained ; 
if we wish to have the latitude itself it will bo given by the 
formula 

latitude = - i* r ' + -//' - &c. 

which, to the degree of approximation adopted, will clearly 
be the same as a. 

lieduction. 

80. We may now consider the term which we had 

neglected (Art. 83) in the expression for the longitude, namely, 

- sin 
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Let N be the position of the node when the moon's 
longitude is 0, M the place of the 
moon, m the place referred to the 
ecliptic. 

Therefore Tw = 




The right-angled spherical triangle NMm gives 

, r tan Nm 

cos Al = XT , 

tan NM ' 



. 1 cosJV tan 2V3f tan 2Vwi 

therefore - iy = .TTT- 

1 + cos N tan NM + 



or, since both JV and NM Nm are small, 



tan'iV NM- Nm . A . 

- , = - . - - approximately ; 
4 rr 



therefore JVJ/-^7//i=^/; 3 sin 2(^ 7)=^ A 3 sin 2(^7), nearly. 

Hence this term, which is called the reduction, is approxi- 
mately the difference between the longitude in the orbit and 
the longitude hi the ecliptic. 



RADIUS VECTOR. 



90. To explain the physical meaning of the terms in the 
value of u. 

We shall, for the explanation, make use of the formula 
which gives the value of u in terms of the true longitude, 
Art. (52). 
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Neglecting the periodical terms, we have for the mean 
value 

tt = a(l-J*'-4wi f > 

The term ^m*, which is a consequence of the disturbing 
effect of the sun, shews that the mean value of the moon's 
radius vector, and therefore the orbit itself, is larger than if 
there were no disturbance. 

Elliptic Inequality. 

91. To explain the effect of the term of the first order ^ 
u = a }1 + e cos (c0 a)}, 



This is the elliptic inequality, and indicates motion in an 
ellipse whose eccentricity is e and longitude of the apso 
a + (l c)0; and the same conclusion is drawn with respect 
to the motion of the apse as in Art. (73). 

Evcction. 

92. To explain the physical meaning of the term 
V mea cos {(2 - 2m - c) 6 - 2/3 + a} . 

This, as in the case of the corresponding term in the longi- 
tude, is best considered in connection with the elliptic 
inequality, and exactly the same results will follow. 

Thus, calling )), O, and a' the true longitudes of the moon, 
sun, and apse (the last calculated on supposition of uniform 
motion), these two terms may be written, 

-a')- 



where E cos S = e -f V me cos 2 (a' o), 

E sin S = i/we sin2 (a - o). 
These are identicaUwith the equations of Art, (77). 



SO LUNAR TIIEOKY. 

Variation. 



93. To explain the cffwt of f/ic term in 1 a cos {(2-2i) 6- 2#}, 

u = a [I -f m* cos [(2 - 2mj 0-2/3}] 
= a[l + 'cod2Q}- )J. 

As far as this term is concerned,* the moon's orbit would 
be an oval having its longest diameter in quadratures and 
least in syzygies. Princlpin, lib. r. prop. 60, cor. 4. 

The ratio of the axes of the oval orbit will be 

. 2 = H !) nearly, m being '0748. 
See lYtncijriit) lib. in. prop. 28. 

Reduction. 

94. The last important periodical term in the value of u is 

73 

cos2(</# 7). 

This term expresses approximately the variation in the 
diifcrence between the values of u in the orbit and in the 
ecliptic. 

For if tt t be the reciprocal of the value of the radius 
vector in the orbit, 

n l = u cos (latitude), 
= =:M ( 1 - 4A "early; 



therefore u w t = i us 9 = ^/; 3 sin 3 ( <;5 7) 

= a& - JaA* cos 2 (g0 - 7) 
= const. \ak* cos 2 (<7^ 7). 

Taking this result in connexion with that of Art. (89) 
wo see that the reduction, so far as periodical terms to the 
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second order are concerned, is simply a geometrical con- 
sequence of the inclination of the orbit ; and that, if we 
measured the longitude and the radius vector along the 
orbit instead of taking their projection on the piano of the 
ecliptic, these periodical terms would not appear. 

95. There are no thcr terms of the second order in the 
value of u. The annual equation, which, in the longitude, 
is of the second order, is only of the third order in the 
radius vector, 

Periodic, time of the Moon. 

96. We have seen, Art. (72), that the periodic time of 
the moon is greater than if there were no disturbing force ; 
but this refers to the mean periodic time established on an 
interval of a great number of years, so that the circular 
functions in the expression are then extremely small com- 
pared with the quantity pt which has uniformly increased. 

When, however, we consider only a few revolutions, these 
terms may not all be neglected. The elliptic inequality and 
the evcction go through their values in about a month,' the 
variation and reduction in about half-a-month ; their effects, 
therefore, on the length of the period can scarcely be con- 
sidered, as they will increase one portion and then decrease 
another of the same mouth. 

But the annual equation takes one year to go through its 
cycle, and, during this time, the moon has described thirteen 
revolutions ; hence, fluctuations may, and, as we shall now 
shew, do take place in the lengths of the sidereal months 
during the year. 

We have, considering only the annual equation, Art. 

pt = + ?>me sin (m0 + -?). 

M 
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Let T be the length of the period: then, when 6 is 
increased by 2-rr, t becomes t+ T; 

therefore p (t+T) = 27T+ 6 + %me! sin (Smir 4- m0 + /3 - ), 
whence ^ ^^ 2<7r + 6 6 ' sin W7r cos OT 4- w# + /3 ) 5 



therefore T mean period -f sinwtTr COS(Q - ), 

where o = uiQ -f /S -f mir = sun's longitude at the beginning 

of the month -f ?^?r 

= sun's longitude at the middle of the month. 
Hence Twill be longest when - =0, 
and shortest when O - = TT ; 

or T will be longest when the sun at the middle of the month 
is in perigee, and shortest when in apogee ; but, at present, 
the sun is in perigee about the M of January, and in apogee 
about the 1st of July; therefore, owing to annual equation^ 
the winter months will be longer than the summer months, 
the difference between a sidereal moiuh in January and July, 
from this cause, being about 20 minutes. 

97. All the inequalities or equations, which our expres- 
sions contain, have thus received a physical interpretation. 
They were the only ones known before Newton had estab- 
lished his theory, but the necessity for such corrections was 
fully recognized, and the values of the coefficients had already 
been pretty accurately determined ; still, with the exception 
of the reduction, which is geometrically necessary, they were 
corrections empirically made, and it was scarcely to be ex- 
pected that any but the larger inequalities, viz. those of the 
first and second orders which we have here discussed, could 
be detected by observation : we find, however, that three 
others, have, since Newton's time, been indicated by obser- 
vation before theory had explained their cause. These are 
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the secular acceleration , discovered by llalley ; an inequality, 
found by Mayer, in the longitude of the moon, and of which 
the longitude of the ascending node is the argument ; and 
finally an inequality discovered by Burg, which has only of 
late years obtained a solution. For a further account of 
these, as also of some other inequalities which theory has 
made known, see Appendix, Arts. (108) to 
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APPENDIX. 

In this chapter will be found collected a few propositions 
intimately connected with the results or the processes of the 
Lunar Theory as explained in the previous pages. Refer- 
ence has been made to some of them in the course of the 
work, and the interest and importance of the others are 
sufficient to justify their introduction here. 

98. The moon is retained in her orbit In/ the force of 
gravity, that is, l>y the same force which acts on bodies at 
the surface of the earth. 

The proof of this is merely a numerical verification ; the 
data Acquired from observation are, 

the space fallen through from rest in 1" by bodies at the 

earth's surface = 16*1 feet, 

the radius of the earth = 4000 miles, 

the periodic time of the moon = 27J days, 

the distance of the moon from the earth's centre=GOx4000 miles. 

The force of the earth's attraction oc r -^ . Therefore, the 

(dist.) 3 * 

space fallen through in 1" at distance of moon by a body 

1 *"! 

moving from rest under the earth's action = - - feet 

-00447 feet. 



MOON'S ORBIT CONCAVE. 8;> 

But the moon in one second describes an angle k ~ -.- *=&>, 
during which the approach to the earth 

= CO x 4000 x 5280 (vers. CD) feet 
00 x 4000 x 5280 x 2?r 2 



= -00418 feet. 

Therefore, the space through which the moon is deflected 
in one second from her straight path, is just the quantity 
through which she would fall towards the earth, supposing 
her to be subject to the earth's attraction ; and we may, 
therefore, conclude that she is retained in her orbit by the 
force of gravity. 

When first Newton, in 16(56, attempted to verify this 
result, he found a difference between the two values equal 
to one-sixth of the less : tho reason of his failure was the 
incorrect measures of the earth, which he made use of in 
his computation; and it was not till about 16 years later 
that he was led to the true result, by using the more 
correct value of the earth's radius obtained by Picart. 

lib. in., prop. 4. 



99. The moon's orbit is everywhere concave to the sun. 

Let $, Ej and M be the centres of the sun, earth, and 
moon. We ought first to apply 
to each body forces equal and 
opposite to those which act on 
the sun in order to bring him to rest. These forces are, 
however, so small that we may neglect them, and we shall 
consider the moon to be moving about the sun fixed, and to 
be disturbed by the earth alone. 
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The forces on M are, therefore, -<.,.,, ........ , ..... in J/$, 

oJ*. 

and ri/:J ............... in ME. 

J j\-L 

This List must be resolved into two, one in MS, the other 
perpendicular to it. 

Therefore, the whole central force on the moon in MS 

7/1 E 



and the proposition will be proved if we shew that this force 
is always positive. JNow, Art. (25), 



m 



. i. . (, /7T.AMJ UTT. O.J/ , 

periodic time or sun = ;I = ,, nearly, 

1 'ttt * in 4 ^ ' 

and moon 

therefore ,, _- 3 j ^ A , lf3 nearly ; 

therefore ^, /3 > T J- u - 

wi' E 
therefore - ,, - - is positive : 

but this is the value of the central force corresponding to 
cosif=-l, and is therefore its least value. Hence the 
central force always tends to the sun, and the path is always 
concave. 

At new moon the force with which the moon tends to the 
sun is, therefore, greater than that with which she tends to 
the earth: the earth being itself H motion in the same 
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direction, and, at that, instant, with greater velocity, will 
easily explain how, notwithstanding this, the moon still 
revolves about it.* 

Central and Tangential Disturbing Forces. 

100. We have hitherto considered the effects of the 

central and tangential disturbing forces in combination ; but 

it will be interesting to determine to which of them the 
several inequalities principally owe their existence. 

(1) To determine the effect of the central disturbing force. 
Make T=(); 

fl v u P 

therefore ' ,, -f u - .-. = 0, 

au~ It u~ 

p 

or substituting for -'. & , from Art. (52), 



2 r 

l* H- w = o | 



1 - :]/t 2 - 1 in 3 + -| nfe cos c0 l - m* cos (2 - 2m) 6, 
-f J wiV cos (2 - 2w - c) O l 



therefore, 

n _ p_l W 

aj 



If we compare this with the value of w found Art. (52), wo 
see that the elliptic inequality and the reduction ar& duo 
to the central or radial force, as also one-half of the varia- 
tion and three-tenths of the evection. 

It would perhaps be proper to separate the absolute cen- 
tral force from the central disturbing force ; the terms due 
to the latter are those which contain m ; therefore, the elliptic 



* .So, the AuUior's Asttroiiowy, p. 23 
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inequality and the reduction are the effects of the former, 
except that in the elliptic inequality the introduction of c, or 
the motion of the apse, is due to the "disturbing force. % 

(2) To determine the rjfect of the tangential disturbing force. 

Let the central disturbing force te zero ; 

P u. 

then p , = Ti(l **) = ) neglecting the inclination, and, 
ft u/ it 

substituting for ISI , . 7/J and for 2( ir ,, +u] L.. a dd their 
// a' av \au J J/iur 

values from Art. (52), the equation in u becomes 

7}J 

. " + u = j I - jj w* cos (2 - 2?//) t + ^ m*e cos (2 - 2m - c) 0J ; 

. f I -f c cos (r0 - a) -f i ?^ s cos f (2 - 2?w ) ^ - 2/3} 

whence ?t = a ^ i i/ rt x /i rt o i 

( -f f ^ we co [(2-2 m - c) - 2/8 + a} . 

We have here the remaining half of the variation and the 
remaining seven-tenths of the evection as the effects of the 
tangential disturbance. Also c = l, or, to, the second order, 
the tangential force has no effect on the motion of the apse. 

101. The separate effects of the central and tangential 
forces in producing the inequalities in the longitude may be 
traced in a similar manner. We shall leave this as an 
exercise for the student. The annual Aquation will be found 
to be due to the central force. The other inequalities will 
be divided as in the radius vector; except as regards the 
variation, of which four-elevenths are due to the central 
and the remaining seven-elevenths to the tangential force. 

To calculate the value of g to the third order. 

102. We must here make use of the results which the 
approximations to the second order have furnished 5 but as 
the value of g is determined by that term of the differential 
equation whose rrgument is yff 7, we need only consider 
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those terms which by their combinations will lead to it with- 
out rising to a higher order than the fourth. 
From Arts. (51), (31), wo obtain 

s = k mifj6 l + *mk sin (2 - 2m - <j) 6^ 

~~ = " 3wA ' l + cos (2 " 



From these we obtain 

Pa S 
AV ' 

T ds 



5 
a -f- s = Oj to the second order ; 

therefore 2 ( * + J I , a 3 f W = 0, to the fourth order. 

Therefore the equation in s becomes, so far as these terms 
are concerned, 

d ^ -f = (- 'X* + iflX* * a 9 ^^' 

assume s = k sin^fl,, 

therefore 4(1- // a ) = - '2 m'/; -f ^ wi** 



To find the value of c to tlw third order. 

103. Proceeding in a similar manner from the equations, 
Arts. (52), (31), we get 



N 
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, = - !' J. {1 + 3 C08 (2 - 8m) *,}, 



~ 3 = - 
we obtain 



{1 - 3e coac0 l ...-*gme cos (2 - 2m - c) 
= a -f ^w^a {3c -f y^wz 



-dm* cos(2~2m) 6 r l ^me cos(2-2w-c)0J 
x sin (2 - 2i) ^ 
= -gm v {sin (2 - 27w) ^ - 2e sin (2 - 2m - c) t - 



n 3 sin (2 - 2w) ^ -f *gme sin (2 - 2m - c)^} ; 
therefore - 8 . ^ =. 

the other term is of the fifth order. 



={!-- 3wi a cos (2 - 2m) ^ -f **m?e cos (2 -2m-c)^ ; 



therefore 2 r/^ + M ) \p--*d6 
and the equation in M becomes 



Assume w = a(l -f e cosc^-f...); 

therefore ae (1 - c 3 ) = (| ? 2 e -f *?w?e] 

whence c *- 1 - m 3 *m\ 
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104. Hence, to the third order of approximation, 

mean motion of apse _ 1 c f mf -f *fm* S + 75m 
mean motion of node "" g 1 "" f i* fom* ~~ 8 3m * 

and since 771 = ^ nearly, we see that the moon's apse 
progredes nearly twice as fast as the node regredes. 

In the case of one of Jupiter's satellites, the periodic 
time round Jupiter is only a few of our days, and the 
periodic time of Jupiter round the sun is 12 of our years, 
therefore m, the ratio of these periods, is very small. 

Hence, the apse of one of Jupiter's satellites progredes 
along Jupiter's ecliptic, with pretty nearly the same velocity 
as the node regredes, assuming these motions to be due to 
the sun's disturbing force; they are, however, principally 
due to the oblateress of the planet. 

Parallactic Inequality. 

105. In carrying on the approximations to a higher order, 
it is found that the expressions for the forces contain terms 
whose argument is the moon's elongation^ i.e. the difference 
of longitude of the sun and moon. These terms appear 
in P and in T of the fourth order, and since 1 TTC, the 
coefficient of 0, is near unity, the terms will become of the 
third order, and therefore of considerable importance in the 
values of u and of t. We shall work out these expressions 
from the earliest steps ; and we select these terms on account 
of the peculiar use which has been made of them to determine 
the sun's parallax, whence they have received the name :>f 
parallactic inequality. 

Let us go back to the expressions for the forces 
Arts. (23), (24). In the first place we may remark that 
the terms we are about to compute must be independent 
of s ; because any term into the composition of which s enters 
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will necessarily have <j in the coefficient of the argument. 
We may therefore at the outset suppose the orbit to coincide 
with the ecliptic. We have then 

SM' = r" - Zr'MO cos (6 - ff] + MG\ 



a . , r 

~ 6 } ~ * 7"~ + * r 



80 M 

We must now substitute these values in 

p Ji'-i-Jf ,/J/O 7^\ ,,/ i i 
P= 



and pick out the terms which have for argument (0 6'). 

T , (MO EG\ xl . . 
In mi 01M 4- OT1 ., there will be a term 
b&/ 

E<F .. .,. * m y E-M , -,. 
m - /t - - cos(0-0')= - ^j/ cos ( fl -^)- 

In Wl ' r> f oj/s o) c o 3 (^-^') wc Sua ll have the terms 



"' - {I cos(0 - 0') + / cos (0- 0') cos2 (6 - 0')}, 

which produce 

3 .- m ^m'r*E- 

cos (0-0') = 33 -^ .. 



and finally V [ Y ,^ a - - .,., 1 sin(0- ff) will produce 
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Therefore the new terms will be 

P Q m'u* EM lf . ,. _ ., ,EM 

- 



in . . .... s < in (fl _ 0- )= _ * m , 8ln ( e_ ff ). 

h'u? b h'u" E+M ^ J * a A+J/ v ; 

; 
Biuce ~ (= - 4 J (F nearly) is of the second order, these terms 

are of the fourth order. 

106. On the principle of the superposition of small 
disturbances, we may compute the effects of these terms by 
themselves as disturbers of the elliptic path. 



AV 
T 



i) 



du f a \ 

-j-Q = ae sin (c& a) ; 

therefore 

-jy-s j = 0, to the fouitii order, 

- -j. u=a. to the second order, 
dv 



Substituting in the differential equation for u, we get 



Assume 



w = a l-f ......... -f A fl ^' - cos{(i-i)5 /3}-r... ; 

1 jL>-f-iJz (t J 
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~ V^ 8 * 

therefore A = - - - - ^ = - 1 g w ; 

1 (1 r/i) 

therefore w=a f"l + ...... - |Jm J~^^ COS{(I-TW) 0-/3J+ ... I . 

107. The corresponding term in the value of 6 will also 
be of the third order, 






/T 
7 , rf^ is of the fourth order and will not rise in t : 
fl'U 

therefore 

* 1 f -' 

39-pL + ...... 



~ pt + ...... ^m -rt .r~ sin {(l i)jp j8} +.... 

Mayer was the first who applied this term to the deter- 
mination of the sun's parallax, by comparing the analytical 
expression of the coefficient with its value as deduced from 

V 

observation. The values of m and of -^=., and therefore of 

M 7 

T? - Kf n' 

_- being pretty accurately known, will be determined, 

that is, the ratio of the sun's parallax to that of the moon ; 
but the moon's parallax is well known ; therefore, also, that 
of the sun can be calculated. The value so obtained by 
Lplace was 8*6"; bui Hansen, by tbo discussion of a great 
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number of modern observations, has found 8*92", a value 
which agrees very closely with recent determinations by other 
methods. 

Secular Acceleration. 

108. Halley, about 1693, found, by the comparison of 
ancient and modern eclipses, that the moon's mean revolution 
is now performed in a shorter time than at the epoch of the 
recorded Chaldean and Babylonian eclipses. This pheno- 
menon, called the secular acceleration of the moon's mean 
motion^ has not yet been fully accounted for. For a long 
time it was altogether unexplained, but in 1787 Laplace 
gave what seemed to be a satisfactory and complete account 
of it. 

The value of p, Art. (54), on which the length of tho 
mean period depends, is found, when the approximation is 
carried to a higher order, to contain the quantity e tho 
excentricity of the earth's orbit. Now, this exccntricity is 
undergoing a slow but continual change from the action of 
the planets, and therefore p, as deduced from observations 
made in different centuries, will have different values. 

The value of p is at present increasing, or the mean 
motion is being accelerated, and it will continue thus to 
increase for a period of immense, but not infinite duration; 
for, as shewn by Lagrange, the actions of the planets on 
the excentricity of the earth's orbit will be ultimately 
reversed, e' will cease to diminish and begin to increase, and 
consequently p will begin to decrease, and the secular 
acceleration will, so far as this cause is concerned, become a 
secular retardation. 

It is worthy of remark that the action of the planets on 
the moon, thus transmitted through the earth's orbit, is more 
considerable than their direct action. 
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109. In his investigation, Laplace treated the excentricity 
e of the earth's orbit as a constant in the differential equa- 
tions, and only considered its variation in the final result. 
This of course greatly simplified the work, but ought to have 
been looked upon only as a first approximation ; yet, strange 
to say, the result so obtained agreed almost perfectly with 
that which had been deduced from - a comparison of ancient 
and modern observations. " Malheureuscment cette iner- 
vcilleusc concordance a laquelle Laplace a du pcut-etre une 
partic dc Fe*clat de sa plus belle dccouvcrte, est sensiblement 
alt^ree par les approximations suivantcs."* 

Some years ago Professor Adams took up the question 
ab initio without any limitation, and introduced the variability 
of e into the differential equations themselves. This strictly 
correct process gave a value of the secular variation only 
about half of Laplace's which had agreed so well with 
observation. 

The variation of excentricity does not therefore account for 
the whole observed change, and the mean motion of the moon 
is affected by some other cause or causes at present unknown. 
Perhaps the remaining acceleration is only apparent, and 
anses from a gradual retardation of the earth's diurnal motion. 
The ciFect would be precisely the sr e, and such a retar- 
dation may be due to the action of the tides. Whatever be 
the cause, the explanation of this unexplained phenomenon is 
a question worthy of the mathematician's most serious efforts. 

Inequalities depending on tie Figure of the Earth. 

110. The earth, not being a perfect sphere, will not at- 
tract as if the whole of its mass were collected at its centre : 
hence, some correction must be introduced to take into 
account this want of sphericity, and some relation must exist 
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between the oblateness and the disturbance it produces. 
Laplace in examining its effect found that it satisfactorily 
explained the introduction of a term in the longitude of the 
moon, which Mayer had discovered by observation, and the 
argument of which i the true longitude of the moon's 
ascending node. 

By a comparison of the observed and theoretical values 
of the coefficient of this term, we may determine the 
oblateness of the earth with as great accuracy as by actual 
measures on the surface. 

111. By pursuing his investigations, with reference to 
the oblateness, in the expression for the moon's latitude, 
Laplace found that it would there give rise to a term in 
which the argument was the true longitude of the moon. 

This term, which was unsuspected before, will also serve 
to determine the earth's oblateness, and the agreement with 
the result of the preceding is almost perfect, giving the 
compression ^^?* which is about a mean between the 
different values obtained by other methods. 

Perturbations due to Venus. 

112. After the expression for the moon's longitude had 
been obtained by theory, it was found that there was still 
a slight deviation between her calculated and observed places, 
and Burg, who discovered it by a discussion of the observa- 
tions of Lahire, Flamsteed, Bradley, and Maskelyne, thought 
it could be represented by an inequality whose period would 
be 184 years and coefficient 15". This was entirely con- 
jectural, and though several attempts were made, it was not 
accounted for by theory. 

* Pont&oulant, Systcme du Mcide, vol. iv. 
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About 1848, Professor Hansen, of Seeberg, in Gotha, 
having commenced a revision of the Lunar Theory,- found 
two terms, which had hitherto been neglected, due to the 
action of Venus. One of them is direct and arises from 
a * remarkable numerical relation between the anomalistic 
1 motions of the moon and the sidereal motions of Venus 

* and the earth ; the other is an indirect effect of an inequality 
'of long period in the motions of Venus and the earth, 
' which was discovered some years ago by the Astronomer 
Royal. 1 * 

The periods of these two inequalities are extremely long, 
one being 273 and the other 239 years, and their coefficients 
are respectively 27'4" and 23*2". 'These are considerable 

* quantities in comparison with some of the inequalities already 

* recognised in the moon's motion, and, when applied, they 
' are found to account for the chief, indeed the only 
'remaining, empirical portion of the moon's motion in 
4 longitude of any consequence ; so Aat their discovery may 

* be considered as a practical completion of the Lunar Theory, 
' at least for the present astronomical age, and as establishing 
*thn entire dominion of the Newtonian Theory and its 

* analytical application over that refractory satellite. 'f 

Motion of tfie Ecliptic. 

118. We have seen, Art. (14), that our plane of reference 
is not a fixed plane, but its change of position is so slow 
that we have been able to neglect it, and it is only when 
tlie approximation is carried to a higher order that the 
necessity arises for taking account of its motion. 



* Report to the Annual General Meeting of the Royal Astronomical Society, 
JFcb. It, 18 W. 

f Address of Sir John Hersohel to the Meeting of the Royal Astronomical 
Socit.y. 
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It has been found to have an angular velocity, about an 
axis in its own plane, of 48" in a century, and the correction 
thus introduced produces in the latitude of the moon a term 
- ceo cos (0 <), 

where o> is the angular velocity of the ecliptic, - the angular 

c 

velocity with which the ascending node of the moon's orbit 
recedes from the instantaneous axis about which the ecliptic 
rotates, </> the longitude of this axis at time I, and 6 the 
longitude of the moon at the same instant. 




Let TAm be the position of the ecliptic at time fy 

A the point about which it is turning, YA = <, 
MN the moon's orbit, M the moon, and Mm a perpen- 
dicular to the ecliptic ; Twt = 6 ; Mm = lat. = 0. 
i the inclination of the orbit, and N the longi- 
tude of the node. 

Let APN'm be the ecliptic after a time St. 
Any point whose longitude is L may be considered as 
moving perpendicularly to the ecliptic with a velocity 
<y sin ( <). 

Hence, the point N will move in the direction PN with 
a velocity o> sin (Nfy. And N 1 will move along PN' with 
a velocity o> sin (N- #) coti ; 

therefore - - = sin (N $) coU*. 

Again, the. point of the ecliptic 90 in advance of N 
will move towards the moon's orbit with a velocity 
sin (90 



100 LUNAR THEORY. 

therefore -r = w cos (N~ <). 

AT 1 d ( N - 4>) l 1, 

Now, cot t, o), and - T - ' = may be consiuciuu v*m- 
777 d* c J 

stant in integrating ; 

therefore SN c& cos (JV 0) cot t, 

&' = cw sin (^V- <), 
and if JV3f = ^, we have 



. . . 

cos sin 't v 
Now, sin /3 = sin i. sin ty ; 
cos /8 . S/3 = cos t . sin >Jr . 81 + sin i. cos - 



but cosi\sin^==cos/3sin(0-J\ r ) and cos^ = cosy8cos(0 N) ; 
therefore S/3 = ow cos (0 0). 

The discovery of this term is due to Professor Hansen; 
its coefficient is extremely small, about 1*5"; but, being of 
a totally different nature from those due to successive 
approximations, it was thought desirable to examine it, and 
the above investigation, which was communicated to me by 
J. C. Adams, Esq.,* will be read with interest on account 
of its elegance. 

With respect to tho foregoing investigation, perhaps the following remarks 

c/iV 
will not be superfluous : The value of * is not the actual velocity of N t 

but its velocity relatively to the position of the node as determined when tho 
motion of the ecliptic is neglected ; its integral is therefore $N the change of 
longitude due to this motion, and in this integration no constant is added, zero 

being taken for the mean value. The periodic forma of both - - and &V shew 
that they oscillate about mean values, the time of a complete oscillation being 

* Now Lowndcan Professor of A&tronomy in the University of Cambridge. 
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that required by sin (^V </>) and cos (*V </) to go through their cycle. This 
relative motion of the node is analogous to that of a particle moving in a straight 
line under the action of a force varying directly as the distance. 

Similar remarks apply to the inclination. 

It must also be borne in mind that the result obtained tJ/3 is not the difference 
between the latitude referred to tho actual ecliptic and that referred to a fixed 
plane, but is the difference between the calculated latitude referred to the actual 
ecliptic on supposition of its being fixed, and the correct latitude referred to the 
same actual ecliptic when, its ration is taken into account. 

We may obtain an approximate value of the coefficient 
ceo by substituting for it , where a is the number of seconds 

through which the ecliptic is deflected in one year = 0'48", 
and n is the number of years in which the node of the moon's 

orbit makes a complete revolution = 18'G ; for then, ~- is 
the angle described by the node in one year; therefore 
- is the ratio of o> : - - -,- -- , supposing < to remain con- 
stant, which is nearly the case ; therefore 

no. 9-3x0-48" ,. ,,* 

CCO = = == 1 4: & . 

27T 3-14 

., 

* This affords the solution of a problem proposed in the Scnafcc-IIou&o 
in the January Examination ^. Id5-. Question *21, Jan. 22. 
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CHAPTER VIII. 

HISTORY OF THE LUNAR PROBLEM BEFORE NEWTON. 

114. The idea which most probably suggested itself to 
the minds of those men who first considered the motion of 
the moon among the stars, was that it described a circle with 
uniform velocity about the earth as a centre ; and this first 
rough result is represented in our equations by neglecting 
all small terms and writing 6 =p, u = a. 

It must, however, have been very soon perceived that 
the actual motion is far from being so simple, and that 
the moon moves with very different velocities at different 
times. 

115. The earliest recorded attempts to take into account 
the irregularites of the moon's motion were made by 
Hipparchus (140 B.C.). lie imagined the moon to move with 
uniform velocity in a circle, of which the earth occupied, 
not the centre, but a point nearer to one side. By a similar 
hypothesis ho had accounted for the irregularities in the 
sun's motion, and his success in this led him to apply it also 
to the moon. 

It is clear that, on this supposition, the moon would seem 
to move faster when nearest the earth or in perigee, and 
slower when in apogee, than at any other points of her 
orbit, and thus an apparent unequal motion would be 
produced. 
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103 



Let BAM be a circle, CA a radius, E a point in A G 
near (7; Cfi, ED two parallel 
lines making an angle a with 
CA. *{ 

Suppose a body Jfto describe J?y - 
this circle uniformly with an ^ 
angular velocity p } the time being 
reckoned from the instant when 
the body was at J7, and the longi- 
tude as seen from E being reckoned 
from the line ED ; 
therefore DEM = 0, BCM = 




EG 
Now - -' ' is a small fraction, and if we represent it by 

we shall have 

lif AO . t -mtr 

sin Jlf= /r ,^sm -4JFJ/ 
O./1/ 

= e sin (0 a) 

/- a), 



l ----- -- -~ 

1 ecos (7^ a) 

this would give ifef, an' 1 then by the formula 

This was called an excentric, and the value of e was called 
the excentricitij) which, for the moon, llipparchus fixed at 
sin 5 1'. 

116. Another method of considering the motion was by 
means of an epicycle, which led to the same result. 

A small circle PJlf, with a radius equal to EC of previous 
figure, has its centre in the circumference of the circle BPD 
(which has the same radius as that of the ex centric), and 
moves round E with the uniform angular velocity p, the 
body M being carried in the circumference of the smaller 



104 



LUNAR THEORY. 




circle, the radius PM remaining parallel to itself, o*, which 

is the same thing, revolving 

from the radius PE with the 

same angular velocity /;, so 

that the angle EPM equals 

PEA. 

Now, when the angle AEP 
equals the angle ACM of the 
former figure, it is easily seen 
that the two triangles EP 
ECM arc equal, and there- 
fore the distance EM and the angle AEM will be the same 
in both, and the two motions will be identical. 

117. The value of e being small, we find, rejecting e*, &c., 

M=c sin (/> a), 
therefore = pt + e sin (pt - a). 

If we reject terms of the second order in our expression 
for the longitude, and make c = 1, we get, Art. (55), 

0=pt + 2e sin (p*-a), 

which will be identical with the above if we suppose the 
cxccntricity of the excentric to be double that of the elliptic 
orbit. 

Ptolemy (A.D. 140) calculated the excentricity of the 
moon's orbit, and found for it the same value as llipparchus, 
viz. 

sin 5 I'^j nearly. 

The excentricity in the elliptic orbit is, we know, about ^ , 
and ^ and ^ will pretty nearly reconcile the two values of 
given above. This shews us, that for a few revolutions the 
inoon may be considered as moving in an excentric, and her 
positions in longitude calculated on this supposition will be 
correct to the first order. 
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Her distances from the earth will not however agree ; for 
the ratio of the calculated greatest and least distances would 



-- - ---- 



, while that of the true ones would be ---- ^ 

l - 



be -- or 
A 



It would, therefore, have required two different excentrics 
to account for the changes in the moon's longitude and in 
her radius vector. Changes in the latter could not, however, 
be easily observed with the rude instruments the ancients 
possessed, and it was very long before this inconsistency was 
detected. 

118. We have said that the moon's longitude, calculated 
on the hypothesis of an exojutric, will be pretty accurate for 
Si Jew revolutions. 

The data requisite for this calculation are, the mean 
angular motion of the moon, the position of the apogee, and 
the magnitude of the cxccntricity. 

But it was known to ilipparchus and to the astronomers of 
his time, that the point of the moon's orbit where she seems 
to move slowest, is constantly changing its position among 
the stars. Now this point is the apogee of Ilipparchus'e ex- 
centric, and he found that he could very conveniently take 
account of this further change by supposing the excentric 
itself to have an angular motion about the earth in the same 
direction as the moon, so as to make a complete revolution 
in about nine years, or about 3 in each revolution.* 

This motion of the apsidal line follows also from our 
expression for the longitude, as shewn in Art. (75). It is 
there, however, connected with an ellipse instead of an 



* On the supposition of an epicycle, this motion of the ap,e couM as 
be represented by supposing the radius which connects the rnoori with the centre 
of the epicycle to have this uniform angular velocity <-f about .> rj in each revolu- 
tion, and dl-o if the -.amc dmclio'i. i 
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cxccntric ; and though the discovery that the elliptic ns the 
true form of the fundamental orbit was not the next in the 
order of time after those of Hipparchus, yet, as all the irre- 
gularities which were discovered in the intervening seventeen 
centuries are common both to Ilippa^chus's excentric and to 
Kepler's ellipse, it will be as well for us to consider at once 
this new form of the orbit, 

Elliptic Form of the Orbit. 

119. We need not dwell on the steps which led to this 
great and important discovery. Kepler, finding that the 
predicted places of the planet Mars, as given by the circular 
theories then in use, did not always agree with the computed 
ones, sought to reconcile these variances by other combina- 
tions of circular orbits, and after a great number of attempts 
and failures, and eight years of patient investigation, he 
found it necessary to discard the cxccntrics and epicycles 
altogether, and to adopt some new supposition. An ellipse 
with the sun in the focus was at last his fortunate hypothesis? 
which waa found to give results in accordance with obser- 
vation ; and this form of the orbit was, with equal success, 
afterwards extended to the moon: biu the departures from 
elliptic motion, due to the disturbing force of the sun, are, 
in the case of the moon, much greater than the disturbances 
of the planet Mars by the other planets. 

In Kcplers's hypothesis, then, the earth is to be considered 
as occupying the focus of an ellipse, in the perimeter of which 
the moon is moving, no longer with either uniform linear or 
angular velocity, but in such a manner that the radius vector 
sweeps over equal areas in equal times. 

This agrees with our investigation of the motion of two 
bodies, Art. (10). 
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Evection. 

120. The hypothesis of an excentric, whose apse line has 
a progressive motion, as conceived by Hipparchus, served to 
calculate with considerable accuracy the circumstances or 
eclipses; and observations of eclipses, requiring no instru- 
ments, were then the only ones which could be made with 
sufficient exactness to test the truth or fallacy of the sup- 
position. 

Ptolemy (A.D. 140) having constructed an instrument, by 
means of which the positions of the moon could be observed 
in other parts of her orbit, found that they sometimes agreed, 
but were more frequently at variance with the calculated 
places; the greatest amount of error always taking place at 
quadrature and vanishing altogether at syzygy. 

What must, however, have been a source of great per- 
plexity to Ptolemy, when he attempted to investigate the 
law of this new irregularity, was to find that it did not 
return in every quadrature, in some quadratures it totally 
disappeared, and in others amounted to 2 39', which was 
its maximum value. 

By dint of careful comparison of observations, he found 
that the value of tins second inequality in quadrature was 
always proportional to that of the first in the same place, 
and was additive or subtractive according as the first was 
so ; and thus, when the first inequality in quadrature was at 
its maximum or 5 1', the second increased it to 7 40', which 
was the case when the apse line happened to be in syzygy 
at the same time.* 

* It would seem as if Hipparchus had felt the necessity for sorao further 
modification of lus first hypothesis, though he was unable to determine it; for 
there is an observation made by him on the moon in the position here specified 
when the error of his tables would be greatest ; and at a time also when she 
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But if the apse line was in quadrature at the same time as 
the moon, the second inequality vanished as well as the first. 

The mean value of the two inequalities combined was 
therefore fixed at 20^'. 

121. To represent this new inequality, which was sub- 
sequently called the 7'7m,Y/o//, Ptolemy Imagined an exeentric 
in the circumference of which the centre of an epicycle moved 
while the moon moved in the circumference of the epicycle. 

The centre of the exeentric and of the epicycle he supposed 
in syzygy at the same time, and both on the same side of 
the earth. 




Thus, if E represent the earth, 

S * sun, 

M moon, 

o the centre of the exeentric EKT in syzygy, 
jR, the centre of the epicycle, would also be in 

syzygy. 

Now conceive r, the centre of the exeentric, to describe a 
email circle about E in a retrograde direction cc', while jR, 
the centre of the epicycle, moves in the opposite direction, 



was in tlic nonagesimaJ, so that any error of longitude, arising from her yet 
uncertain parallax, would be avoided. Ptolemy, who recorda the observation, 
employs it to calculate the cvcction, and obtains a result agreeing with that 
f)i his own observations. (Sec Delambre. Ast t Aaciennc.) 
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in such a manner that each of the angles S'Ec, S'ER may 
be equal to the synodical motion of the moon, that is, her 
mean angular motion from the sun ; SES' being the motion 
of the sun in the same time. 

Now we have seen, Art. (110), that the first inequality 
was accounted for by supposing the epicycle EM to move 
into the position rw, r and li being at the same distance 
from i T , and rm parallel to 11 Al* the first inequality being 
the angle rEm. But when the centre of the epicycle is at 
72', and 11' M' is parallel to rtn y the inequality becomes R'EM' } 
and we have a second correction or inequality iuEM'. 

122. That this hypothesis will account for the phenomena 
observed by Ptolemy, will be readily understood. 

At syzygics, whether conjunction or opposition, the centres 
of the exccntric and epicycle are in one line with the earth 
and on the same side of it ; the points r and li' coincide, as 
also m and M'. Hence mEM' = 0. 

At quadratures (figs. 1 and 2) c and R are in a straight 
line on opposite sides of the earth, and therefore It' arid r at 

1%. 2. 





their furthest distance. If, however, M' and m be at the 
same time in this line, or, in other words, if the apse line 



* For simplicity we leave out of considerate n the motion of the apse. 
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be in quadratures (fig. 1), the angle inEM' will still be zero, 
or there will be no error in the longitude. But, if the apse 
lino is in syzygy (fig. 2), the angle mEM' attains its greatest 
value.* 

Ptolemy, as we have said, found tins greatest value to be 
2 39', the angle mEr being then 5 1'. 

123. Copernicus (A.D. 1313), having seen that Ptolemy's 
hypothesis gave distances totally at variance with the obser- 
vations on the changes of apparent diameter,! made another 
and a simpler one which accounted equally well for the in- 
equality in longtitude,]and was at the same time more correct 
in its representation of the distances. 

Let E be the earth, OD an epicycle whose centre G de- 
scribes the circle C' GO" about E with the moon's mean 
angular velocity. 




Let CO, a radius of this epicycle, be parallel to the apse 



* If Ptolemy had used the hypothesis of an excentric instead of an epicycle 
for the first inequality of the moon, an epic\cle would have represented the 
second inequality more simply than his method did. Dr. Whewell's History of the 
Inductive *SVtVwc'tv, vol. I., p. 230. 

f See Delauibre, Ast. Mafcrnf, vol. I., p. 116. Whewell's History of Inductive 
Sciences, vol. I., p. 395. 
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line EA, and about as centre let a second small epicycle 
be described, the radii CO and OMbc so taken that 

CO- OM . pQ . .C0+0J/ . . A , 

n~v sin 5 l ) an< * 7^7^ = sin 7 40 . 

OH \jTj 

The radius OM must now be made to revolve from the 
radius OC twice as rapidly as EC moves from ES, so that 
the angle COM may be always double of the angle CE8. 

From this construction, it follows that in syzygies the 
angle CE8 being or 180, the angle COM is or 360; 
and therefore C and M are at their nearest distances, as in 
the positions C' and C"' in the figure. Then CM= CO - OM, 
and the angle CEM will range between O e and 5 1', the 
greatest value being attained when the apse line is in quad- 
rature. 

When the moon is in quadrature CES=QQ or 270, 
and, therefore, G Y #J/=180 or 540 and C and M are at 
their greatest distance apart, as in the position C" ; then, 
CM CO + OM, and ihe angle CEM will range between 
and 7 40', the former value when the apse line is itself 
in quadrature, and the latter when it is in syzygy. 

124. Thus the results attained by Ptolemy's construction 
are, as far as the longitudes at syzygics and quadratures are 
concerned, as well represented by that of Copernicus; and 
the variations in the distances of the moon will be far more 
exact, the least apparent diameter being 28' 45" and the 
greatest 37' 33" ; whereas, Ptolemy's would make the greatest 
diameter 1.* 

The values which modern observations give vary between 
28' 48" and 33' 32". 



* Delambre, Ast. Mvderne. 
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125. It will not now be difficult to shew that the *intn>- 
duction of this small epicycle corresponds with that of the 
term l ine sin {(2 e 2m - r) pt - 2/8 -f aj in our value of 9. 

For, referring to the preceding figure, we have 

OEM= sin OEM = ~ mi OME 



= ~ sin (2 . 8 KG - A EM) 

OM . , , , . , , v 

= / y sm i- (moon s mean long. - sun s long.) 

(moon's true. long. - long, of apse)}, 

and OEM being a small angle whose maximum is 1 19J-', 
we may write moon's mean longitude instead of the true 
in the argument, and also EC for OE* therefore, 

OEM . * sin {2 (moon's mean longitude sun's longitude) 
(moon's mean longitude longitude of apse)} 



= 4770" sin {(2 - 2wi - c)jp< - 2/3 + a}. 

The value of the coefficient is from modern observations 
found to be 4589-61". 

120. In Art. (77), we have considered the effect of this 
second inequality in another light, not simply as a small 
quantity additional to the first or elliptic inequality, but as 
forming a part of this first; and therefore, modifying and 
constantly altering the excentricity and the uniform pro- 
gression of the apse line. 

Boulliaud (A. D. 1645), by whom the term Evection was 
first applied to the 'ccuntl inequality, seems to hint at some- 
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thing of this kind in the rather obscure explanations of his 
lunar hypothesis, which, never having been accepted, it would 
be useless to give an account of it.* 

In Ptolemy's theory, Art. (121), the evection was the 
result of an apparent increase of the first lunar epicycle 
caused by its approaching the earth at quadratures ; but, in 
this second method, it is the result of an actual change in 
the elements of the elliptic orbit. 

D'Arzachel, an Arabian astronomer, who observed in 
Spain about the year 1080, seems to have discovered the 
unequal motion of the apsides, but his discovery must have 
been lost sight of, for Ilorrocks, about 1640, re-discovered it 

* in consequence of his attentive observations of the lunar 

* diameter : he found that when the distance of the sun from 

* the moon's apogee was about 45 or 225, the apogee was 
c more advanced by 25 than when that distance was about 
c 135 or 315. The apsides, therefore, of the moon's orbit 
4 were sometimes progressive and sometimes regressive, and 

* required an equation of 12 3()', sometimes additive. to their 

* mean place and sometimes subtractive from it.'f 

Ilorrocks also made the excentricity variable between the 
limits -06686 and '04362, * 

The combination of these two suppositions was a means 
of avoiding the introduction of Ptolemy's excentric or the 
second epicycle of Copernicus : their joint effect constitutes 
the evection. 



* Aprds avoir etabli lea mouyements et les epoquea de la lime, Boulliaud revient 
& 1'explication de Vcvection ou de la seconde inegalite. Si sa theorie u'a pas fait 
fortune, le nom du moms cat rcste. ' En merne temps que la lune avancc sur son 

* cone autour de la terre, tout le gynteme de la lune tut d'/tlace ; la terre (important la 
'lune, rejette loin d'elle 1'apogee, et rapproche d'autant le perigee; mais cette 

* (jyection a des bornes fixees.' Delambre, Hist, de F Ast. Mod., torn. 11. p. 157. 

f Small's Astronomical Discoveries of Krpler, p. 307. 

Q 
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Variation. 

127. After the discovery of the evection by Ptolemy, a 
period of fourteen centuries elapsed before any further ad- 
dition was made to our knowledge of the moon's motions. 
Hipparchus's hypothesis was found sufficient for eclipses, and 
when corrected by Ptolemy's discovery, the agreement be- 
tween the calculated and observed places was found to 
extend also to quadratures; any slight discrepancy being 
attributed to errors of observation or to the imperfection 
of instruments. 

But when Tycho Brah6 (A. D. 1580) with supei'ior instru- 
ments extended the range of his observations to all inter- 
mediate points, he found that another inequality manifested 
itself. Having computed the places of the moon for differ- 
ent parts of her orbit and compared them with observation, 
he perceived that she was always in advance of her com- 
puted place from syzygy to quadrature, and behind it from 
quadrature to syzygy ; the maximum of this variation taking 
place in the octants, that is, in the points equally distant 
from syzygy and quadrature. The moon's velocity therefore, 
so *ar as this inequality was concerned, was greatest at new 
and full moon, and least at the first and third quarters.* 



* 'It appears that Mohammed- Aboul-Wefa-al-Bouedjani, an Arabian astro- 
nomer of the tenth century, who resided at Cairo, and observed at Bagdad in 975, 
1 discovered a third inequality of the moon, in addition to the two expounded by 

* Ptolemy, the equation of the centre and the evection. This third inequality, the 

* variation, is usually supposed to have been discovered by Tycho Brahe, six 
'centures later In an almagest of Aboul-Wefo, a part of which exists in 

* the Royal Library at Paris, after describing the two inequalities of the moon, 
'he has a Section IX., "Of the third anomaly of the moon called Muhaznt or 

'jProjMMiMiV 1 But tliis discovery of Aboul-Wefa appears to have excited 

' no notice among his contemporaries and followers ; at feast it had been long 

* quite forgotten, when Tycho Brahe re-discovered the same lunar inequality.' 
WhewelTs Hist, of Inductive F'ienccs, vol. i. p. 243. 
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Tycbo fixed the maximum of this inequality at 40' 30". 
The value which results from modern observations is 
39' 30". 

128. We have already two epicycles, or one epicycle and 
an excentric, to explain the first two inequalities: by the 
introduction of another epicycle or excentric, the variation 
also might have been brought into the system ; but Tycho 
adopted a different method:* like Ptolemy, he employed 
an excentric for the evection, but for the first or elliptic 
inequality be employed a couple of epicycles, and this 
complicated combination, which it is needless further to 
describe, represented the change of distance better than 
Ptolemy's. 

To introduce the variation, he imagined the centre of the 
larger epicycle to librate backwards and forwards on the ex- 
centric, to an extent of 40^' on each side of its mean position ; 
this mean place itself advancing uniformly along the excentric 
with the moon's mean motion in anomaly ; and the libration 
was so adjusted, that the rnoon was in her mean place at 
syzygy and quadrature, and at her furthest distance from it 
in the octants, the period of a complete libration being half 
a sy nodical revolution. 

Annual Equation. 

129. Tycho Brah<$ was also the discoverer of the fourth 
inequality, called the annual equation. This was connected 
with the anomalistic motion of the sun, and did not, like the 
previous inequalities, depend on the position of the moon : i 
her orbit. 



* For a full description of Tycho's hypothesis, see Delambre, Hist, de ? Ast. 
Mod. torn. I. p. 102, and An Account of the Astronomical Discoveries of Kepler 
by Robert Small, p. 130. 
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Having calculated the position of the moon corresponding 
to any given time, he found that the observed place was 
behind her computed one while the sun moved from perigee 
to apogee, and before it in the other half year. 

Tycho did not state this distinctly, but he made a cor- 
rection which, though wrong in quantity and applied in an 
indirect manner, shewed that he had ^een the necessity and 
understood the law of this inequality. 

He did not try to represent it by any new excentric or 
epicycle, but he increased by (8m. 13s.) sinew's anomaly] 
the time which had served to calculate the moon's place;* 
thus assuming that the true place, after that interval, would 
agree with the calculated one. Now, as the moon moves 
through 4' 30" in 8m. 13s., it is clear that adding (8m. 13s.) 
sin (surfs anomaly) to the time is the sa^ie thing as sub- 
tracting (4' 30") sin (suns anomaly) from the calculated lon- 
gitude, which was therefore the correction virtually intro- 
duced by Tycho.f Modern observations shew the coefficient 
to be 11' 0". 

We have seen, Art. (82), how this inequality may be 
inferred from our equations. 

Seduction. 

130. The next inequality in longitude which we have to 
consider is not an inequality in the same sense as the fore- 
going r that is, it does not arise from any irregularity in the 
motion of the moon herself in her orbit, but simply because 
that orbit is not in the same plane as that in which the 



* That is, the equation of time which he used for the moon differed by thai 
quantity from that used for the sun. 

f liorrocks (lti"ii>) made the correction in the same manner as Tycho, but sc 
increased it that the corresponding coefficient was 11' 51" instead of 4' 30" 
Flainstced was the first to r pply the correction to the longtitude instead of the time 
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longitudes are reckoned, so that even a regular motion iii 
the one would be necessarily irregular when referred to the 
other. Thus if NMn be the moon's 
ecliptic, and if M the moon be re- 
ferred to the ecliptic by the great 
circle Mm perpendicular to it, then 
MN and mN are (f 90, 180, 270, 
and 360 simultaneously, but they 
differ for all intermediate values : the difference between 
them is called the reduction. 

The difference between the longitude of the node and that 
of the moon in her orbit being known, that is the side NM of 
the right-angled spherical triangle NMm^ and also the angle 
N the inclination of the' two orbits, the side Nm may be 
calculated by th<j rules of spherical trigonometry, and the 
difference between it and NM^ applied with a proper sign 
to the longitude in the orbit, gives the longitude in the 
ecliptic. ^ % 

Tycho was the first to make a table of the reduction 
instead of calculating the spherical triangle. Jlis formula 
was 

reduction = tan 2 ^/ sin2L ^ tan 4 ^/ sin4i, 

where /is the inclination of the orbit and L the longitude of 
the moon diminished by that of the node. 

The first term corresponds with the term %k 2 sin2 (gpt 7) 
of the expression for 0. 

Latitude of the Moon. 

131. That the moon's orbit is inclined to the ecliptic wa& 
known to the earliest astronomers, from the non-recurrence 
of eclipses at every new and full moon ; and it was also 
known, sinces' the eclipses did not always take place in the 
same parts of the heavens, that the line of nodes represented 
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by Nn, in the preceding figure, has a retrograde motion on 
the ecliptic, AT moving towards T. 

Hipparchus fixed the inclination of the moon's orbit to 
the eclipcic at 5, which value he obtained by observing the 
greatest distance at which she passes to the north or south 
of some star known to be in or very near the ecliptic, as for 
instance the bright star Regulus ; and by comparing the re- 
corded eclipses from the times of the Chaldean astronomers 
down to his own, he found that the line of nodes goes round 
the ecliptic in a retrograde direction in about 18$ years. 

This result is indicated in our expression for the value 
of the latitude by the term k sin (yd 7), as we have shewn 
Art. (85). 

132. Tycho Brand further discovered that the inclination 
of the lunar orbit to the ecliptic was not a constant quantity 
of 5 as llipparchus had supposed, but that it had a mean 
value of 5 8', and ranged through 9' 30" on each side of 
this, the least inclination 4 58^' occurring when the node 
was in quadrature, and the greatest 5 17^' being attained 
when the node was in syzygy.* 

lie also found that the retrograde motion of the node was 
not uniform : the mean and the true positions agreed very 

* Ebn Jounis, an Arabian astronomer (died A.D. 1008), whose works were 
translated about 50 years since by Mons. Sedillot, states that the inclination of the 
moon's orHt had been often observed by Aboul-Hassan-Aly-ben-Ainajour about 
the year 918, and that the results he had obtained were generally greater than 
the 5 of Hipparchus, but that they varied considerably. 

Ebn Jounis adds, however, that he himself had observed the inclination 
Bcveial times and found it 5 3', which leads us to infer that he always observed 
in similar circurastanceg, for otherwise a variation of nearly 23' could scarcely 
have escaped him. Sec Delambre, Hist, de V Ast. du Moyen Aye, p. 139. 

The mean value of the inclination is 5 8' 55-40", the extreme values are 
4 57' 22" and 5 20' 6". 

The mean daily motion of the line of nodes is 3' 10-64", or one revolution in 
6793-29 days, or 18 y. 218 d. 21 h. 22m. 46 s. 
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well when they were in syzygy or quadrature, but they 
were 1 46' apart in the octants. 

By referring to Art. (87) we shall see that these correc- 
tions, introduced by Tycho Brahe, correspond to the second 
term of our expression for s. 

Since Hipparchus could observe the moon with accuracy 
only in the eclipses, at which time the node is in or near' 
syzygy, we see why he was unable to detect the want of 
uniformity in the motion of the node. 

133. To represent these changes in the position of the 
moon's orbit, Tycho made the following hypothesis. 

Let ENF be the ecliptic, K its pole, BA C a small circle, 




having also K for pole and at a distance from it equal 
to 5 8'. Then, if we suppose A the pole of the moon's 
orbit to move uniformly in the small circle and in the 
direction BAG } the node N^ which is at 90 from both A 
and K, will retrograde uniformly on the ecliptic, and the 
inclination of the two orbits will be constant and equal to A K. 
But instead of supposing the pole of the moon's orbit to 
be at -4, let a small circle abed be described with A as polo 
and a radius of 9' 30" ; and suppose the pole of the moon's 
orbit to describe this small circle with double the velocity 
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of the node in its synodical revolution which is accomplished 
in about 346 days, in such a manner that when the node is 
in quadrature the pole may be at a, the nearest point to K, 
and at c the most distant point when the node comes to 
syzygy, at d in the first and third ojtants, and at l> in the 
second and fourth, so as to describe the small circle in about 
173 days, the centre A of the small circle retrograding mean- 
while with its uniform motion. 

By this method of representing the motion, we see that 

when node is in quadrature) . ,, / f A~rt=58' 9l'=458j-', 

1 ^tionoi tlie-^ 

.................... y*ygy J orbit is Uc= 

while at the octants it has its mean value 



Again, with respect to the motion of the node, since N is 
the pole of KaAc, it follows that when in syzygy and quad- 
rature, the node occupies its mean place ; hi the first and 
third octants, the pole being at d, the node is before its 
mean place by the angle dKA = (9' 3D") coscc5 8' 1 46', 
nearly, and it is as much behind in the second and fourth 
octau cS. 

So that the whole motion of the ncue, and the correction 
which Tycho had discovered, were properly represented by this 
hypothesis, which is exactly similar to that which Copernicus 
had imagined to explain the precession of the Equinox. 



SELECTION OF EXAMINATION QUESTIONS 

FROM COLLEGE AND SENATE-HOUSE PAPERS AND FROM THE 
MODERATORSHIP \ND FELLOWSHIP EXAMINATIONS AT 
TRINITY COLLEGE, DUBLIN. 

1. Define the plane of the ecliptic and prove that, as 
seen from the earth, sin X=; T ^^ (TiT sinV nearly; X being 
the latitude of the sun and A/ that of the moon. 

2. Obtain the differential equation of the moon's radius 
vector. 

3. What is the principle on which the successive approxi- 
mations to the moon's path are obtained ? 

Compare the coefficients of the principal term in the cen- 
tral disturbing force with the principal term in the central force. 

4. If in solving the* equation 

fu =a _P _ T c//<_ 
dOP + U ~1tu* h'u 8 d0 

we take as our first approximation a a {l+c cos(0-a)}, 
we obtain at our higher approximations terms containing in 
their coefficient. Kow is this detect avoided by taking 
u = a (1 + e cos(c# a)} as our first approximation? 

Upon what principle do we approximate to the value of c ? 

T 

5. Shew that in the expansion of ^ -.. there will be a 

nu 

term - ^ wiV sin {(2 -2m- 2r) 0-20 + 2aj which w^uld 
rise to the second order in the longitude. 

Shew further that the term with the same argument does 
not appear actually in the longitude assuming 

)+...+ V ecos {(2- 2m- c) 0-20 + a.} 
V me 2 cos {(2 - 2m -2c) 0-20 + 2a}+.. .]. 
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G. The longitude of the moon contains a term of the form 

A sin |(2 - 2m - c] pt -2/3-1- a} 

where A is of the second order; find how it must have 
arisen a^d determine A. 

7. Investigate the following expression for the moon's 
longitude as far as the second order, supposing the orbits 
of the moon relative to the earth, and of the earth relative 
to the sun, to be originally circles in the plane of the ecliptic, 

=pt + y in* sin {(2 - 2m) pt - 2/3}. 

8. Suppose that whatever constants are involved in the 
arguments 0, </>, ^, &c. in the equation 

L = A -f B sin 6 -f G sin $ -f I) sin ^r -f &c. 
are perfectly known, and that observation gives the numerical 
values of a considerable number of longitudes L^ L^ L^ vKrc., 
corresponding to the known angles tf,,^,^,, &c., 0,,0,,^.,, &c., 
a , $ a , tjr !t , VY.C., Ac., shew how the constants A, />, 6 f , &c. 
may be nunierically dcterminod in the two distinct cases 
where the observations at our disposal arc and are not 
unlimited in number. 

9. KXJ lain how this process fails and may be modified 
when two of the terms in question Ji sin 0, G sin <f> are nearly 
synchronous in their periods. 

10. Shew that the evection in longitude, viz. 

tfnw sin |(2 - 2m - c) pt - 2/3 + ctj, 

may be represented as the joint effect of certain periodic 
changes in the exccntricity of the lunar orbit and in the 
mean longitude of its apse. 

11. Assuming the usual notation of the Lunar Theory, 
explain the physical meaning of the following equations : 

s = A sin(# 7), 

s~k sin (00 7), 

s = A- sin (ffd - 7) + g ink sin {(2 - 2m - g) B - 2/8 + 7 J . 

12. Explain the effect of the term i\i cos{(2-2w) 0-2$)} 
in the moon's radiui vector; find the number of days in 
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the period of the resulting inequality and the ratio of tho 
axes of the oval orbit. 

13. Investigate the effect of the annual equation on the 
length of the lunar month. 

Assuming v^='075, v "017, sin l.T.'SOWaS.'Ui), and tho mean 
period of a sidereal revolution of the moon = 27d. 7h. -Mm. IrJs., 
find the difference between a winter and a summer month. 

14. Shew that the moon's orbit in space is always concave 
towards the sun. 

Assuming that the earth describes a circle about the sun, 
and the moon a circle about the earth and in the same piano 
with the earth's orbit, compare the curvatures of the moon's 
orbit in space at her perihelion and aphelion. 

15. Prove that if we go to the third order of approxi- 
mation, the motion of the moon's apse in one revolution 
of the moon equulj '-^nf (1 V'w) 300. 

P T 

10. The expressions for ., , and ,. ,, contain terms of 
1 h it, h u 

the form 3 A cos (#--#') inJ. A sin(0 9'} respectively; com- 

pute the terms resulting from them in the values of a and 0. 

t 

17. Copernicus represented the evection by an epicycle 
superposed on the epicycle which represented the elliptic 
inequality. From the construction given by him sjicw 
geometrically that the evection vanishes when the moon's 
mean elongation from the sun is half her true anomaly, 
and, in general, varies nearly as the sine of twice tho 
difference of these angles. 

18. The secular change of inclination of tho actual ecliptic 
to a fixed ecliptic being a" annually, shew that this will 
give rise to an inequality in the moon's latitude whose typo 






is t y~a" cos<, nearly; n being the number of years in 

7T 

which the moon's nodes make a complete revolution, and 
< the difference of longitude of the moon and of the ascending 
node of the actual ecliptic. 

\v. MKK \i.ri; AM* ^"V, I'liivrcus rvMiu:iJ>iiK 
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